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The rate of growth of an organism may be described as a function of both its mass and of the 
time. In particular, if the rate of growth is described as a function of mass alone, one derives the 
so-called principle of “equifinality,”’ according to which the final mass achieved is independent 
of any interruption of the growth process. Since this phenomenon is observed in many in- 
stances, the assumption that the rate of growth is a function of mass alone has some theoretical 
justification. On the other hand, it is obvious that in many instances irreversible effects can be 
achieved by interrupting the growth process, thus indicating that the rate of growth may 
depend on variables other than mass. The consequences of some of these alternative assump- 
tions are examined. A growth factor G is postulated, which is supposed to decay with time, and 
whose concentration determines the over-all rate of growth per unit mass. The rate of decay of 
G is calculated from the postulated growth equation under two alternative assumptions, (1) 
_ that Gis a constituent of all the cells of the organism and (2) that G is produced separately and 
distributed through the organism. It is shown that in each case the concentration of G satisfies 
a simple differential equation. While in normal growth the results deduced from these hypothe- 
ses are indistinguishable from those deduced from the assumption that the rate of growth 
depends on mass alone, widely different results are obtained in the case of interrupted growth. 
In particular, the rate at which growth is resumed after interruption and the final mass 
achieved depend on the time and on the duration of interruption. It appears, therefore, that 
various assumptions about the nature of the “aging” of the organism lead to various quantita- 
tive relations between growth and the times of interruption and resumption of growth. 


It has been noted that, during the growth of an organism, the masses of 
several organs remain proportional to constant powers of the total mass 
(and therefore to powers of each other). Such growth, called “allometric,” 
is described by the system 


m;=a;M"*, (a;, a; > 0) (1) 


where m; is the mass of the 7th organ, M the total mass, a; and a; con- 
stants characteristic of the ith organ. On the other hand, the time course 
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of the total mass seems in many cases well described by an empirical 
equation proposed by L. von Bertalanffy (1951) 


ON = Me — kM (ight 0: 2 al (2) 


Bertalanfty ascribes the following theoretical interpretation to equation 
(2). The positive term on the right represents the rate of anabolism (the 
building up of tissues), and the negative term the rate of catabolism (the 
breaking down of tissues). The fact that in many organisms u = 3 gives 
further theoretical significance to this interpretation if one supposes that 
the rate of anabolism is proportional to the total surface and that of 
catabolism to the total mass of the organism. This assumption seems espe- 
cially reasonable for spherical bacteria, for which the value » = 3 typi- 
cally holds. Moreover, in insects (in the larval stage) the value w = 1 is 
typical, which is again reasonable in view of the direct air profusion of 
tissues of these small organisms. The same value, u = 1, holds for many 
rod-like bacillae, as one would expect of organisms growing principally in 
one direction. Finally, equation (2) implies that the rate of growth of the 
total mass depends only on the mass and not on time explicitly. This con- 
clusion is borne out by some arrested growth experiments, where, after the 
growth of an organism is allowed to resume after interruption, it is re- 
sumed at a rate approximately equal to the rate at which it had been 
interrupted, so that the organism “‘catches up” and achieves its normal 
final mass. Some of these experiments are cited by L. B. Mendel (1923), 
Bertalanffy (loc. cit.), and by R. W. Gerard (1940) and discussed as indica- 
tive of the possibility that arrest of growth is also an arrest of “aging.” 

The actual reversibility of the growth process, as in planaria (Gerard, 
loc. cit.), supports a suspicion that at least some organisms are actually 
“ageless.” However, the question of the reversibility or irreversibility of 
biological processes is an extremely involved one. Irreversible arrests of 
growth can certainly be achieved, not necessarily through such obviously 
irreversible procedures as the removal of vital glands in higher organisms, 
but also by processes whose irreversible character is not so pronounced. 
The ancient practice of foot binding among the Chinese is one example 
which comes to mind. Retarded growth based on malnutrition in humans 
is also notoriously irreversible if carried beyond the period of normal mat- 
uration. Hence it is not unreasonable to postulate the working of some 
“aging mechanism” at least in higher animals. We undertake to investi- 
gate the theoretical consequences of some such postulates, assuming that 
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the normal growth curve can be empirically described by equations (1) 
anid (2). 
Equations (1) and (2) can be combined into a system of differential 
equations 


ae a;m,f 

(3) 
dM 
a rpm 


The quantity f = dM/Mdt is, of course, the rate of growth per unit mass 
of the total organism by definition. Various hypotheses can now be made 
about the nature of f. Whatever hypotheses are made, however, the sys- 
tem (3) states that the rate of growth per unit mass of each organ is pro- 
portional to the same quantity f, the constant of proportionality being 
characteristic of each organ. Call these constants, a;, the allometric con- 
stants, taking as unity the allometric constant of the whole organism. 

When the allometric equations are put into the form given by (3), atten- 
tion naturally centers on f, which appears in all equations. If fis a function 
of the total mass alone, i.e., 


f= AM —F, (4) 


as equation (2) implies, growth cannot depend on age explicitly and will 
be resumed after an interruption exactly as it had stopped. If, however, 
we confine our attention to normal growth, as given by the solution of (2) 
for M as a function of time, then the form of f is by no means prescribed. 
In particular, f can be taken to be an explicit function of ¢, such that the 
solution for M(¢) in (3) will coincide with that in (2). To be sure, this 
assumption about f will not lead to the phenomenon of equifinality and 
thus will contradict the experimental evidence cited. We will nevertheless 
pursue the consequences of this assumption for their intrinsic interest and 
for the questions to which they give rise. 

We seek a function f(#), such that M(é), as obtained from the solution 
of (3), will be identical with M(é), as obtained from a solution of (2). 
Separating the variables in (2), we obtain the integral 


M dM 
ale hie — kM’ >) 
The quadrature can be obtained in closed form and solved for M, namely, 


, : 1/(1—n) 
= |—_— —k(1—p) t 
M F Ke »| We, 66) 
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where 


K =4—-M-«(0) : es 


Formula (6) is a generalization of Bertalanffy’s solution for the special case, 
where » = 2. The case » = 1 leads to exponential growth and will be 
considered separately. In what follows y will be assumed to be less than 
unity. 

Since f, considered as a function of ¢, must be identically equal to the 
function of M as given by (4), we find that 


poh{t— Kew] 8 (8) 
where K is given by (7). 

We now turn our attention to f. We note that the units of f are those of 
k, namely, [7]. This is, of course, a consequence of the definition of f as 
the time rate of change of mass per unit mass. If we wish to interpret f as a 
concentration (expressed as mass per unit mass) of a hypothetical growth 
substance G, f must be dimensionless. We accordingly imagine f multiplied 
by a constant denoting the time unit, and thus a unit of f becomes the con- 
centration of G which will produce a unit increase of mass per unit mass 
per unit time. We note that f is monotone decreasing, and this decrease 
must be attributed to a ‘‘cause.” The flow of time alone is not sufficient in 
any complete formulation of a physical situation to account for any 
change. The simplest hypothesis would be to assume that the rate of 


change of f is a function of f alone. Indeed, if we differentiate (8) with 
respect to /, we find 


dj - esl bike es 
dP ary eae fn (9) 
ET Rete | 


while combining terms in (8), we find that 


i Ree eta 
eae ae (10) 
Rk 
h 
ftk= 
4 Keto! (11) 


Thus we have proved 


Theorem 1. Let growth be determined by the concentration { of a growth 
substance G, such that the time rate of change of { follows a course dependent 
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on f alone. Then, if the growth coincides with the solution of equation (2), f 
satisfies the differential equation 


- 
Gea Uwe). (12) 


We have noted that this interpretation of the system (3) is contradicted 
by the phenomenon of delayed growth, because if f decays “automati- 
cally” in accordance with equation (12), while the mass remains constant, 
then the resumption of growth at a later time should start at a rate smaller 
than the rate observed at interruption. If, however, we suppose that the 
suspension of growth suspends also the decay of f [for example, if f is “used 
up” in accordance with (12) in the growth process only], then this conclu- 
sion does not follow. On the other hand, this latter supposition is strictly 
ad hoc. In fact, it is logically equivalent to admitting that f depends on 
mass alone. 

Still, an argument can be made along those lines. The point is that 
suspension of growth can be made in a variety of ways. In the experiments 
cited, growth was suspended by diet deficiencies. It may be that the very 
deficiencies which inhibited growth also inhibited the decay of f. It is cer- 
tainly possible to design methods of ‘“‘growth suspension” which will give 
different results. To mention a very crude method, suppose some ‘‘non- 
essential” but non-regenerating portions of the animal, like the extremi- 
ties, are amputated during a certain period in such a way that the total 
mass remains constant even as the growth of the remaining portions of the 
body proceeds. It is certainly unlikely that the animal will “catch up” and 
achieve normal mass, compensating for the loss of the amputed parts by 
increasing the remaining mass. The question is whether a process can be 
found which would achieve similar results, while the animal appeared 
“normal” in other respects. Such a process would have to “disconnect” G 
(if it exists), thus preventing growth and yet allowing the supposed “ag- 
ing” process, as reflected in the decay of f, to continue. If such a process 
were found, let us see how the interrupted growth curve would look; in 
particular, how the final mass would depend on the times of interruption 
and resumption of growth. 

In thus proceeding, we may make either of two assumptions: (1) we 
may assume that equation (12) represents the decay of f independently of 
the increase of mass; or (2) we may assume that equation (12) includes the 
decrease of f (which is a concentration) due to the increase of total mass. 
Assumption (1) would be in order if G were a constituent part of every cell 
and were reproduced as cells reproduce. In this case the local decrease of 
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G (i.e., the decrease of f) would have nothing to do with the increase of the 
total size of the organism. Assumption (2), on the other hand, would be in 
order if G were thought of as being produced by some organ and diffused 
throughout the organism. Under this assumption, even if the total mass of 
G remained constant, while the mass of the organism increased, f would 
decrease, because G would have to be distributed over a greater mass. At 
any rate, if f continues to change while M is kept constant for a period 
At, then after growth is resumed, we should have instead of dM/dt = 
M (df(t) 
dM 


ay MO FAD . (13) 


We shall deduce the consequences of assumption (1) first. 

In what follows the following notation will be used: 

M(t) : mass at time ¢ under normal conditions of growth. 

M*(t) : mass at time ¢ under conditions of interrupted growth. 

f*(é) : concentration of G under conditions of interrupted growth. 

t, : time of interruption of growth. 

¢, : time of resumption of growth. 
We will also use subscripts instead of arguments of functions, e.g., M1 for 
M(h), M= for M*(), etc. 

To determine the final mass under normal conditions, we set dM/dt = 0 


in (2) and obtain 
h\i- 
Mo= (=) ; (14) 


Now in the interval f; < ¢ < & if f continues to decay according to (12), 
then 


f Kke*Q-+)¢, 
Noe . 
ut K e-*(—n) ty (15) 
k 
On the other hand, 
Mo E Were eHow 4] a (16) 
2 kb . 
From (4) and (14), we get 
Ma]! 
pt mise 
fo Ka (17) 
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However, the designation of the time origin is arbitrary. Let us designate 
tz as the time origin. Corresponding to equation (17), we will then have 


Kke*(-#)t, 
ja fe AABN 
: Soave Ke*(-»t, M; (18) 


whereupon, using (16) and re-arranging, we obtain 


h 


k = Mo (19) 
p—Ketows a rey cane 
hth 
=—-|>-— —k (ln) t, 
eee a hs | ae ae 
aeons ie eee Tie 
poke k (1—p) ty 2 
Vi— 
us =(5) Fm Me (21) 


We have proved 
Theorem 2. Let growth be governed by G, whose concentration obeys (12). 
Then if growth is interrupted between t; and tx, the ratio of the final mass to 
the normal final mass equals the ratio of the normal mass at t; to that at te. 
Similarly the rates of growth can be compared at # and t}. We have 


h 
iscaye 9 etamwy |S — Kew] 


ahd & Sitek a ae nesses (22) 
*7 
M (41) Ai e—F(1—n) t, & — K e-*0-») | 
For normal growth 
1h 
aly (ile = ea Te l= 
ett [eM (d |, oe 


so that (22) reduces to 


M*’ (i) _ [Mo *— Ma "| eee (24) 
MG) (MM 


It is of interest to note that if » = 2, M!“ = M”* is roughly the measure 
of the linear dimensions of the organism. Designating this average linear 
size by L, we have for interrupted growth 


M*’ (tz) a (Leo ~ Ly) Li (25) 
M’ (h;) (Lo — Ly) Le 
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In the growth curve reproduced by Bertalanffy (Joc. cit.), Mi/Mz is 
approximately 4, while M/M. is approximately unity. Thus by the re- 
sult of Theorem 2, the hypothesis herein proposed seems to be dramati- 
cally refuted. Still the question remains whether by some other method of 
growth interruption or in some other organisms results more in accordance 
with equation (21) could not be obtained. 

We now pass to assumption (2), in which account must be taken of the 
increase of total mass as contributing to the decrease of f. Here evidently 
the decay of f during the interval of constant mass will not be according to 
equation (12) but more slow, possibly even negative. We seek the modified 
equation of decay. Let F = Mf be the total mass of G. Then 


dF df dM 


ap ea) tae oa 
Substituting (3) and (12) into (26), we obtain 
dF 
“Tp el aia) al 
(2) 


=M(uf?—fk+ufk). 


If growth is interrupted, M remains constant, and so during interrupted 


growth we have df'/dt = Mdf/dt. Hence during the interval of constant 
mass 

df 

Fiimtie o iek F (4S5¢Sh). (28) 

It is interesting to note that whereas under assumption (1), f’(é) was 

always negative, under assumption (2), F’(é), and therefore f’(¢) during 
the interval of constant mass, may be negative, zero, or positive, depend- 
ing on when the growth is interrupted. Let us investigate the sign of f’(é). 
Since f > 0 for all values of f, it is sufficient to examine the sign of uf, — 
k(1 — pw). By (4), we have f; = hM4-! — k. Therefore 


ph— kp) = what — b= hale — 2) | 
— (29) 
-ie(urr-M). | 


Noting that » — 1 < 0, we see that the expression (29) is negative, zero, 
or positive according to whether M}~™ is greater than, equal to, or less than 
uM". Therefore there is a value of ¢, namely, when M1“ = pM ot, SuGgh 
that if growth is interrupted at just that time, f will remain constant 
throughout the period of interrupted growth [since f’(¢) is a function of f 
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alone]. This situation will be indistinguishable from the one in which 
M'(t) is a function of M alone (as in Bertalanffy’s hypothesis). If, how- 
ever, growth is interrupted at a later time, when M}“ > pM", f will 
continue decaying during the interval of interruption, and growth will be 
resumed at a slower rate, leading to a smaller final mass. If growth is inter- 
rupted at 4; when My“ < pM", f will start accumulating. This situation 
presents a difficulty, however, because in this case f can become infinite in 
finite time. This must be interpreted as a failure of the model. We shall 
therefore restrict ourselves to the case where the resumption of growth is 
not delayed so long that 47>“ becomes greater than M1}™, for it is in that 
case that f(¢) must pass through a singularity. 

The integration of equation (28) is straightforward, leading, under the 
initial conditions on fi, to the expression for f* as a function of ¢, namely, 


Kk2(1—) e#a-we 
(1—p) —Kke*0-44 Khe 


r=; (30) 


The ratio of the resumed to the interrupted rate of growth is, as above, 
fe /fi, ie., after proper substitutions, in terms of the normal masses 


M*'(t) (1—y) Mr *(Mo"— Mp *) 
Mt) (Mi pM) (Me = Ms) 


(31) 


If Mi > uM", then a forteriorti My * > pM" for tg < ~, and the 
ratio (31) is always positive. It is also easy to verify that the ratio is less 
than unity and tends to unity as M{“ approaches uMz*. Since we have 
confined our considerations to the case where Mj" > wMy™, (31) is posi- 
tive even if Mi“ < uM", but in this latter case the ratio is greater than 
one and tends to unity from above as My“ approaches pM". 

It remains to calculate the ratio of the final masses, W3/M... This can 
be done by considering two organisms starting with initial values fs and 
fi respectively of f and with identical masses M > = M, att = 0, and by 
comparing their final masses. From (17) we obtain 


Mine 
k 


(fotk) =Ma", (32) 


whereas, substituting f; and f2 for fo and M7. >» = My, for M,, we obtain 


ME fac h 
MA Hire’ tate 


(33) 
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After elaborate substitutions, we obtain the ratio in terms of the normal 
masses at ¢; and f, and in terms of pw, namely, 


Me (1—p) —MS (Ms I) 
Mine 1 ; (34) 


(oo) 1 —p( — 
Again, when w = 3, we can express (34) in terms of linear dimensions: 


Ms'* (1-4) —Le (2-4) 


Mo Ls (33) 


In the experiment cited by Bertalanfly mentioned above, Mi, M2, and 
M., are approximately 80,160, and 180 grams respectively, while y is 
assumed to be 2. These values make the ratio (34) about $, which is cer- 
tainly not far from unity and well within the accuracy of the experiment 
(since variations of mass due to factors other than growth, such as water 
content or fat reserve, make higher accuracy improbable). Equation (31) 
could not be tested, because it was impossible to estimate from the data 
the rates of growth at the critical times with any reasonable accuracy. In 
order to test the hypothesis herein proposed, it is necessary to have data 
from which rates of growth could be inferred reasonably accurately, and 
moreover it is necessary to be able to compare the results of interrupting 
growth at various times for various intervals. Fairly accurate estimates of 
growth rates could be made from Chart 4, p. 128, in Mendel (op. cit.). 
Unfortunately, this chart pictures a comparison between normal and re- 
tarded (not interrupted) growth, and so equation (31) could not be ap- 
plied. It is noteworthy, however, that the several resumption rates appear 
to be very nearly equal to the normal rates corresponding to the same 
masses. Again Bertalanffy’s hypothesis seems to be corroborated at least 
for the methods employed. Nevertheless the theoretical investigation of- 
fered here should serve to keep open the question of whether in organisms 
other than rats and/or by methods other than dietary deprivation in the 
interruption of growth other growth patterns cannot be obtained. 

For completeness, we will discuss the case where u = 1. Here we must 
have h > k for positive growth. We so assume and accordingly rename 
h — k by h. Now normal growth is exponential. Moreover, f is a constant, 
and thus our assumption (1) is indistinguishable from the assumption that 
Ff depends on mass alone, which is now vacuously true. Assumption (2), 
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however, leads to a different result. Starting again from equation (26), we 
note that now f’(¢) = 0, so that 

BEd Regs a pel Z 

Cee spelicn is (,StSh). (36) 
Now f is certain to accumulate during the time of interrupted growth. To 
exclude the case where f becomes infinite, we must have ft: < + fp, 
Imposing this restriction and proceeding exactly as before, we obtain in 
this special case for the ratio of resumed to interrupted rates of growth 


MP (Gy. 1 
Mh) ath PSA) ba 
or, in terms of masses alone, 
Meo ort 
Mth) 5, Ma (38) 
itor i, 


The final masses cannot be compared, because both become infinite in this 
case. 

It is realized that the considerations herein proposed are far too simple 
and naive to account for the vastly complicated growth processes, espe- 
cially of higher organisms. Yet the frequently observed regularities of 
many growth curves indicate that perhaps some such simplified scheme 
may be a proper starting place for a mathematical theory of growth. At 
any rate, we believe that the assumptions made by Bertalanffy go a step 
beyond the purely empirical descriptions of growth curves and attempts 
to fit them with empirical equations. In particular, our hypotheses lead to 
the consideration of the possibility of actually accelerating growth in cer- 
tain instances by arresting it and to quantitative relations between the 
characteristics of the growth curve and the times at which growth is ar- 
rested and allowed to resume. Thus, the hypotheses invite further and 
more refined experimental procedures designed to support or to refute 
them. Moreover, a method is indicated for formulating any number of 
alternative hypotheses in which the time rate of growth depends on any 
properly selected and manipulable variables. The technique of arresting 
growth is then tantamount to holding some of these variables constant, 
and thus changing the postulated quantitative relations among the others. 

Finally, we note that the allometric growth equations for the several 
organs lead to an analogous bifurcation of hypotheses concerning the na- 
ture of f. We may, for example, assume that the concentration of G is the 
same in all organs, and that the allometric constant a; represents an ‘‘ac- 
tivity coefficient” of G with regard to its effect on the growth of the 7th 
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organ. Or, on the other hand, we can assume that the allocation of G to 
the several organs differs, so that a,f represents the actual concentration 
of G in the ith organ. The two hypotheses again lead to different results, 
equivalent in the case of normal growth, distinct in the case of interrupted 
growth under assumption (2), because the differential growth rates have 
differential effects on the concentrations. To show this, let us write, in- 
stead of equations (3), the time course of growth for the 7th organ as 


ES = Mf i. (39) 


Then our theory applies equally to the 7th organ under the following 
transformations: 

[= yd ata 

k=kia;- (40) 


ee 
Equation (12) remains invariant under this transformation except for 
the factor a;: 


2 hin) a go hae Os GG alae (41) 


Thus under assumption (1), the allometric constants are preserved after 
the interruption of growth. However, equation (28) is not invariant. We 
have, in fact, 


fluf—k(—w)] =a; fil (usa, +1—a:) fi—ks (1-1, (42) 
so that we may have simultaneously 
pple <0 
(uegel—o.)if; — Rifle) =e — BA — a) 0 (43) 


if a; # 1. Therefore it is possible to arrest growth at a time when f’(#) < 0 
while f;(¢) > 0 (or vice versa), thus distorting the proportions of the or- 
ganism. Such distortion would be a consequence of making both assump- 
tion (2) and interpreting the allometric constants as “allocating” f to the 
various organs. It would not be a consequence of either assumption (1) or 
of an interpretation of the allometric constants as “activity coefficients.” 
Finally our theory can be modified by supposing that the relative rate 
of catabolism (the breakdown of tissues) per unit time per unit mass is 
constant (%), so that f, instead of being given by (4), is given by 


h 
f=hM(t) = ; 
i K ek) t (44) 
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Proceeding as before, we now obtain, instead of (12), 


d 
Ja 6-2). (45) 


Since for positive growth we must have f > k, df/dt is negative fort < t, 
as in equation (12). If we proceed under assumption (1), we find that 
Theorem 2 holds precisely, as is intuitively evident, since changing f by an 
additive constant is tantamount to shifting the time origin. Under as- 
sumption (2), however, we are led instead of to (28) to 


ak 2 
“r= uf), (4th) (46) 


which is always positive for f > &, so that under this assumption, an inter- 
ruption of growth always leads to a greater rate of resumed growth and 
consequently to greater final mass. 

One thus has a choice of combining various hypotheses about the nature 
of a hypothetical factor associated with growth to obtain rather widely 
different theoretical results for the case of interrupted growth. 


I wish to thank Dr. R. W. Gerard and Dr. L. von Bertalanffy for read- 
ing and critically discussing this paper. 
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Simple theoretical models are proposed for the study of the interdependence between car- 
diac contraction, arterial pressure, and capillary drainage. The relation between pressure and 
flow is derived for a model of branching distensible tubes taking into account the finite pulse 
wave velocity. Equations are derived both for the case where the pulse wave is non-distorted 
and for the case where the wave is damped and distorted to a limited extent. Following the 
model of J. W. Remington and W. F. Hamilton (1947), the former case is applied to the 
larger arteries. Expressions are developed for the stroke volume, cardiac ejection, and systolic 
arterial storage in both the steady and non-steady states. Expressions for the percentage dis- 
crepancy involved in the computation of these quantities from a single tube model as con- 
trasted with a multi-branched model are derived. For typical cases these discrepancies are 
small and thus credence is lent to the further use of the simpler single tube model which re- 
quires fewer independent parameters. It is also shown that the formulae for stroke volume 
and arterial storage are only slightly sensitive to changes in pulse wave velocities, and that 
for some purposes it would seem permissible to assume an infinite velocity. 

The problem of capillary drainage is discussed, and the consequences of equations developed 
for the case of a distorted wave are shown to compare favorably with published experimental 
data. An approximate boundary condition for capillary drainage is derived. 

Finally, A. V. Hill’s velocity load equation for muscle is used to obtain a first approximation 
for the velocity of cardiac contraction in terms of the initial arterial pressure, the heart radius, 
and the parameters of the heart musculature. It is shown how methods developed for stroke 
volume determination from the pressure contour may be used to estimate the heart and ‘‘air 
chamber’’ parameters. Use of these parameters and those obtained by other independent 
measurements permits the principle variables to be determined numerically. 


Introduction. The blood circulation in larger animals has evolved pri- 
marily because it serves to reduce the distances that would otherwise have 
to be covered by passive diffusion. It performs this active transport by 
means of an energy source which pumps blood through the arterial tree 
and into the capillary beds. The variable most easily measured is the 
pressure; however, from the above point of view, the variable of most in- 
terest is the amount of material transported per unit time. 

In this paper, we shall be interested in the interaction of events taking 
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place in the heart, arterial tree, and in the capillary beds, focusing our 
attention on the arterial pressure only in so far as it is related to the flow 
of material. 

A variety of information has been obtained on isolated muscle fibers, 
on the contracting ventricle, on the arterial tree, and on the flow of blood 
through capillaries. In the intact animal, the variables specifying these 
systems are not all independent, but are related, in principle, through 
equations of constraint. This paper is motivated by a desire to elucidate 
and find approximations to some of these equations of constraint. 

For example, in the latter portion of the paper an attempt is made to 
apply certain ideas about skeletal muscle to the heart. Even though the 
similarity between cardiac and skeletal muscle has been emphasized by 
many authors (e.g., Frank, 1895; Starling, 1918), there have been rela- 
tively few mathematical attempts to apply skeletal muscle concepts to 
the circulation. If the heart is made up of muscle fibers, a knowledge of 
the initial length of these fibers (which is determined by the diastolic vol- 
ume of the heart) should be sufficient to determine the time course of 
contraction, and hence the ejection velocity of the blood from the ven- 
tricle, provided the load-resisting movement of the fibers is known. The 
determination of the load complicates the problem considerably. More 
specifically, this load is a function of the interventricular pressure which 
is dependent upon the aortic pressure. However, the aortic pressure is 
determined primarily by its own static characteristics and the cardiac 
ejection. Thus we are led in a natural manner to the problem of the deter- 
mination of the cardiac ejection and the stroke volume when the arterial 
pressure contour is known. Among other things, the solution depends on 
the rate at which fluid leaves the arterial tree, i.e., the rate that fluid 
enters the arterioles; hence we must also concern ourselves with the prob- 
lem: Given the pressure at the beginning of the arterioles, what is the flow 
into the capillary bed. The flow out of the capillary bed is constant over 
any single heart cycle, and hence does not present an immediate problem. 

In the outline above, the problem has been broken up in a manner 
which is experimentally realistic. It is to be expected that detailed study 
of any of the above categories would reveal problems that are of interest 
for their own sake, and independent of the proposed synthesis. Some of 
these will be indicated or dwelt upon in this report. However, since the 
primary motivation is synthesis, it is expedient to strive for simplicity, 
and for models whose parameters are capable of evaluation. In many in- 


stances the errors attendant upon proposed simplifications will be esti- 
mated. 
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Since this is a first approximation, we shall assume that the contra cting 
left ventricle is roughly spherical with an approximate radius r and that 
its volume V, is proportional to the cube of its radius. The walls of the 
ventricle are made up of layers of muscle fibers oriented in different direc- 
tions. In general, we may expect to find fibers pulling in every direction. 
However, we shall focus our attention solely on the average tangential 
forces and velocities exerted by these fibers. In other words, we assume 
that upon contraction the size of the ventricle diminishes, but the shape 
remains invariant (the average forces and velocities in directions other 
than tangential are considered to be negligible). 

The assumption of a spherical ventricle has appeared frequently in the- 
oretical discussions (cf. Rushmer eé al., 1951, 1952; Burch et al., 1952). 
Because of the high pressures the ventricle must contract against, one 
might expect its walls to “round up.” More specifically, Rushmer et al. 
(1951, 1953) suggest that the right ventricle contracts as a sphere, but the 
left ventricle is primarily cylindrical and that contraction is brought about 
essentially by a reduction in diameter supplemented by a relatively small 
reduction in length. Ignoring this small reduction of length, the equations 
for a cylindrical chamber will be similar in many respects to those of a 
sphere. We shall proceed with the discussion in terms of a sphere and only 
briefly mention from time to time the analogous results for a cylinder. 

Turning now to the velocity of contraction, O. Horwitz (1951) found 
the rate of approach of two markings placed (at varied and arbitrary posi- 
tions) on the surface of the ventricle remains practically constant through- 
out the contractile phase of the cycle. In terms of our model, this obser- 
vation suggests that the average tangential velocity is constant, or, since 
the average length of the contractile element will be proportional to the 
circumference, which in turn is proportional to the radius, we may put 


dr 
eee (1) 
where is a positive constant, and time is denoted by ¢. 

Equation (1) is implied by the Horwitz experiments. Perhaps the 
plausibility of these results is indicated by the following argument. While 
carrying a constant load, muscle contracts at a constant velocity, and the 
change from rest to constant velocity occurs very rapidly in a twitch. 
Further, in a normally beating heart, the signal to contract occurs when 
the interventricular pressure is below the aortic pressure and the valves 
are not open. Consequently, the contraction is isometric until the inter- 
ventricular pressure is large enough to force open the aortic valves. It is 
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likely that by the end of the isometric contraction period most of the 
transients involved in going from rest to activity have died out. During 
the major portion of contraction, the pressure rises, but on the other 
hand the radius falls, and since the tension is roughly proportional to the 
product of the pressure times the radius, the tension or load may remain 
relatively constant, implying a constant velocity of contraction. Accord- 
ing to R. W. Ramsey (1944), during the late stages of an isotonic muscle 
contraction the velocity begins to diminish. But, during the latest stages 
of cardiac contraction both the radius and pressure decline, and thus the 
load is reduced so that it is possible that the heart may practically finish 
its contraction at the same constant velocity. 
Integration of (1) yields 
r=1ro—X (2) 


where 7p is the resting value of r. We also have, for the ventricular volume, 


Veaeiet (3) 
where A is a positive constant, which for an exact sphere is equal to 47. 
If we define our time scale such that t = 0 at the beginning of each beat 
(more specifically ¢ = 0 at the moment the aortic valves open), then the 
rate of fluid outflow from the heart U will be given by 


dV, 
dt © 


U=— 


However, this is applicable only if the valves are open. If the valves close 
at time ¢ = T, (duration of systole) and the length of each cycle is de- 
noted by 7, then for T, < ¢t < T no fluid flows from the heart to the ar- 
terial tree, i.e., from equations (1), (2), and (3) we have 
U=AX(ro— AL)?2, URS eG 
re } (4) 
= 04 Tare 


For a cylinder, the analogous expressions would be 


US AUN (rp N10), Sars; 

=0, to i Dae 

where 7’ is the radius of the cylinder, —\ = dr’/dt, and A’ ~ 2m times the 
length of the cylinder. 


Let us now turn our attention to the question of how a given cardiac 
ejection effects the time course of the arterial pressure, returning in a 
later section to the dependence of \ on heart and arterial characteristics, 
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The relation between pressure and flow in a network of arteries. Suppose 
a pressure contour is recorded from a network of elastic tubes such as that 
found in the circulation. Then it would be reasonable to ask at what rate 
must fluid have been injected into the tubes to cause the given contour. 
Conversely, if the time course of injection is known, we should be able to 
determine the pressure contour it will give rise to. 

The classical approach to this chronic problem in cardiovascular physi- 
ology has been to consider a long distensible tube which has no longitudi- 
nal curvature and no branches. Starting with the Navier-Stokes equa- 
tions, simplifications have been sought which would allow these equations 
to be solved. Drastic simplifications are necessary because, “In case of a 
distensible tube like a blood vessel, the size and shape of the boundary are 
themselves determined by the distribution of velocities and pressures 
which in their turn are determined by the boundary conditions. Problems 
of that nature . . . lead to complicated, usually non-linear integro-differ- 
ential equations or more general equations in functional derivatives. 
Exact solutions of such equations may well be considered as hopeless at 
the present time” (Rashevsky, 1945). On top of these difficulties, H. Bran- 
son (1945) has mentioned complications in viscosity and compressibility 
due to the presence of blood corpuscles. It is very probable that even 
approximate solutions of these equations give insight into local condi- 
tions, such as the dependence of the pulse velocity on the elastic proper- 
ties of the artery (Jacobs, 1953). However, the very nature of the Navier- 
Stokes equations requires a precise point-to-point knowledge of the com- 
plete geometry of the system; so that even if the exact solutions were 
forthcoming it may be reasonable for our purposes to question the sig- 
nificance of replacing the vastly complex architecture of the arterial tree 
by a single tube of infinite length. Even if it is feasible to lump the major 
arterial branches into an equivalent tube, the length of the tube is not long 
when compared with the pressue pulse wave length, and hence the as- 
sumption of a long tube implies a hydraulic impedance matching between 
the tube and the capillary beds, when in reality the pressure in the capil- 
laries is non-pulsatile. In other words, aside from gaining insight into 
localized conditions, we may ask whether the difficulties and complica- 
tions involved in seeking an exact solution are warranted when there is 
no assurance that the model is sufficiently representative of our actual 
system. In this connection, it is also worth while to keep in mind the 
degree of accuracy which can be expected of the interesting experimen- 
tally measured variables in cardiovascular physiology. 

In view of the foregoing, it seems reasonable, as an alternative ap- 
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proach, to attempt to find some approximation which is simple both in 
principle and in practice and which does not rely so heavily on the exact 
specification of intimate geometrical details (which are unknown). It 
would be desirable for this approximation to incorporate as much of the 
known. information as possible. Yet, on the other hand, it may well be 
meaningless to insist on extreme degrees of precision. 

A significant advance in this direction was made by Otto Frank (1899), 
who introduced an “‘air chamber” analogy. In this model all of the arteries 
are treated collectively and likened to a chamber with the following prop- 
erties: (1) pressure waves within the chamber are transmitted with an in- 
finite velocity, and (2) the chamber leaks at a rate which varies linearly 
with the pressure (analogue of capillary drainage). The condition of con- 
tinuity of mass flow applied to the chamber results in the desired relation, 


1 dP IP 
€ ange 8) 


where P is the pressure in the chamber, U is the rate of fluid flow into the 
chamber, W and ¢ are constants, and 1/e is a measure of the distensibility 
of the chamber, i.e., if V is the volume of the chamber, « = dP/dV. (In 
the original Frank theory c = 0. The introduction of a nonzero ¢ is due 
to Gomez, 1937.) This model has been further elaborated and modified 
by many workers (Apéria, 1940; Karreman, 1954). 

Let us denote e/W by 7; then if (4) is inserted into (5) the integral from 
0 tot < T, of the resultant equation will be 


es 
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The order of magnitude of y as given by A. Apéria (1940) is such that 
for ¢ < T, it is possible to ignore terms of higher order than two in the 
expansion of e~” without serious difficulty (e.g., for c = 60, T = 1, yT, 


is of the order of .3). Expanding and retaining only terms less than the 
third degree in ¢ gives 


P(t) =Pot+ [eAdro—y(Po— ) ] t 


' (7) 
in [2 cAdri—% (Py 0) +eAr’r| p 
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Under physiological conditions, blood initially (i.e., at ¢ = 0) flows intu 
the arterial tree faster than it flows out. This implies 
Po eee 


fa 


Adri > 


or, since y = ¢/W, 
eAdro > y (Po— 6); (8) 
hence the linear term in (7) is positive. Multiplying both sides of (8) by 


7/2 yields 
Gann (Ha-o 


and consequently the coefficient of the quadratic term in (7) is negative. 
Thus in this case the systolic portion of the pressure curve is given by the 
parabola (7) which rises to a maximum and whose coefficients are parame- 
ters of the heart and other parts of the circulation. Similar remarks apply 
to the analogous equation for the cylindrical case which is given by 


P(t) =Pot [eA’N ro —y (Po c)] t—4 [yeA’ A’ ro ted’ 2 — 2 (Po — | i? , 


During diastole, T, < ¢ < T and U = 0. Hence introducing the sym- 
bol P, = P(T.), and integrating (5) from ¢ = T, tot < T we have 


P(t) = 6+ (P;— Cc) e~v(t— Ts) : Wn SERIE « (9) 


According to Yamaguchi (cf. Apéria, 1940) equation (9) represents a good 
approximation to the diastolic pressure curve. It would seem reasonable 
to use equations (7) and (9) as approximations to at least the gross quali- 
tative features of the circulation. However, until the methods of deter- 
mining € and y are specified, relatively little can be said about the possi- 
bility of making quantitative predictions from them. We shail return to 
this problem in a later section where we will develop a method for the 
determination of the parameters which specifies that the pressure con- 
tours (7) and (9) correspond to those contours which yield plausible 
values for the stroke volume. 

For the present, we may note that there has been considerable disagree- 
ment in the literature on the method of determining ¢ and y (Apéria, loc. 
cit.), and that the most frequent objection to the “air-chamber theory” 
has been the assumption of an infinite pulse velocity. Apart from any 
physical analogue, (5) may be viewed purely as a phenomenological equa- 
tion. With U = 0, it merely expresses the frequently observed situation 
where, when flow into the arterial system is interrupted, the pressure de- 
cays exponentially. When U # 0 it assumes that superimposed on this 
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tendency for the pressure to decay will be a tendency for the pressure to 
rise at a rate proportional to the rate of fluid flow into the system. 

An important attempt to introduce the finite velocity of wave propaga- 
tion has been recently made by Remington and Hamilton (1945, 1947), 
although their assumption and results have not been stated in analytical 
form. Briefly, these authors proposed a model which divided the arterial 
system into from four to six segments, each segment having distinct pres- 
sure-volume characteristics, and pressure wave velocities. These quanti- 
ties were carefully measured as a function of diastolic pressure in dogs and 
humans, and the results published in tabulated form (Remington and 
Hamilton, loc. cit.). The major assumptions of their model were: (1) the 
pressure wave is transmitted undistorted throughout the arterial bed, 
(2) capillary drainage occurs at a rate proportional to the pressure at the 
origin of the capillary bed, less 20 mm. Hg, and (3) in a steady state the 
net volume gain of the arterial system at the end of systole must equal the 
amount of fluid that drains out of the capillary beds during diastole, which 
follows whenever U = 0 throughout diastole. After recording the ar- 
terial pressure in the aorta as a function of time, they attempted to deter- 
mine the distribution of fluid ejected by the heart at intervals of 10 m. 
sec., by using the above assumptions and referring to the tables of trans- 
mission times and distensibilities. This enabled them to reconstruct the 
cardiac ejection curve and thus determine such things as stroke volume 
and cardiac work. A method for the determination of stroke volume based 
on the above principles, but differing in practical details, has been recently 
advanced by H. R. Warner et al. (1950). The limits of validity of these 
methods have been the subjects of several recent papers (Huggins ef al., 
1952; Longino et al., 1952; Opdyke, 1952, 1953; Remington et al., 1949). 

In the next few sections, we shall develop analytical expressions relat- 
ing the pressure and flow in a network of arteries. We shall study the in- 
fluence of branching and finite pulse wave velocities on the calculation of 
stroke volumes. The problem of a limited amount of wave distortion will 
be touched upon, and the problem of arteriolar drainage dealt with at 
length. 

The equation of continuity. Consider a system of tubular bifurcations, 
such that each branch either drains into another branch or terminates 
in a capillary bed. 

Let there be a total of W branches, and let the ith branch be divided into 
NV; segments (see Fig. 1). The variables in the jth segment of the ith 
branch will be denoted by the use of the subscript ij. Thus its pressure 
will be given by P,;, its volume by V;;. 
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Let each segment be characterized by the parameter E£,;, which is de- 
fined as 
adP;; 


Ei; oe Van 


(10) 
We shall refer to 1/£;; as the distensibility of the ij segment. In general 
E;; will be a function of V;;; however, we shall consider the particular case 
where £,; remains constant over the range of pressure variations encoun- 
tered during one cycle. That this actually occurs is borne out by the data 
of Remington and Hamilton, both on rings and whole segments of ar- 
teries. Physically it would seem preferable to define an elastic coefficient 
by Vi;(dP;;/dV;;); however, it has been found empirically (Remington 
and Hamilton, 1945) that Z£;; remains practically constant even when the 


Fricure 1. Schematic diagram of the zth branch 


tensions of the musculature imbedded in the arterial wall are altered. It 
is common practice to use equation (10) under dynamic conditions, and 
we shall do likewise, although it should be noted that we are ignoring the 
contribution of fluid accelerations to the pressure. 
If fluid enters the 77 segment at a rate U;;, then it will leave at a rate 
U ij11; hence the rate of change of volume of the segment will be given by 
OV Ui — Us sts. (11) 
Multiplying (11) by £;; yields 
CPi dP ;; aViga. 


[er din eee eee) (12) 
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For any fixed i expression (13) represents V; equations in V;-+1 variables 
of the form U;;. Eliminating Ui2, Uis, . .. Uiw; by successive substitutions, 
we arrive at 

_ yh 1gePs 
— 4t EA dt 
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The term Uj represents the rate at which fluid enters the branch, and 
U iy, 41 the rate at which fluid leaves. If the ith branch drains into a capil- 
lary bed, then U;y,41 is the rate of capillary drainage, and we shall denote 
it as v;. On the other hand, if the ith branch is not a terminal one, then 
it must drain into at least one other branch. If, for example, it drains si- 
multaneously into the &th and the mth branches, 


U; N;t1 a Ua U mi . 


Expressions for Uj: and U1 are given by (14), with 1 = k and 2 = m re- 
spectively. Hence, starting with any arbitrary integer 7 and making suc- 
cessive substitutions for Ua and U;,,,, we arrive at 


1. dP 
Eo dl + 9 (15) 


i 


coz 


where U = Uy denotes the rate that fluid is injected into the system. 
Note that equation (15) has been derived primarily on the basis of the 
conservation of matter. We shall refer to it as the equation of continuity. 

Non-distorted waves. Let us now consider the case, similar to the Rem- 
ington-Hamilton model, of a non-distorted, non-damped pressure wave in 
the larger arteries. While this is, of course, not exactly true, these authors 
nevertheless have been able to make reasonable predictions on this basis 
in many instances. 

Characterizing each segment 77 by an arrival time 7;;, such that the 
pressure wave requires 7,;; seconds.to travel from the origin to that seg- 
ment, then the assumption of a non-distorted pressure wave implies that 


Py=Py (t— 743) =P G—T5) (16) 


where P(#) is the pressure contour at the origin. 
Using (16) in (15)fyields 


N 
~ id 
u= > Edin (tet (17) 


If the number of segments of each branch becomes infinite while the 
length of each segment approaches zero, the summation over the index 
j in equation (17) becomes an integral. To carry out this procedure, let 
1; denote the length of the ith branch, 1/Z; its distensibility (i.e., E; is the 
slope of the pressure volume curve for the ith branch), x,; the distance 
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of the jth segment ith branch, from the origin of the branch, and Az;; the 
length of that segment. We then have 


El; 
AX ;; ; 


Byg= (18) 


If 5; 1s the transmission time to the origin of the ith branch, and c; is the 
transmission velocity through the branch, then 


(19) 
Inserting these values into (17) yields 
ee d 
Sr ee ee i i ee sae Pains 
er a b= 8) any + a. (20) 
Upon passing to the limit, (20) becomes 
“1 id Mig 
The integral in (21) may be evaluated by making the substitutions 
vi 
oy aN ee dah | 
(22) 
GPAS yy OP obs 2 ak | 
dik UNE idl ede: 
and, since ¢ is held constant during the integration 
dx=— c,d &;. (235) 


Thus (21) becomes 


Scola dee 24 
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Let b; + (J;/c:) = e:, the arrival time at the end of the ith branch. Then 
expression (24) for U becomes 


Petey, — Pi — ¢€,) | oe. (25) 


Defining 


AP, =P (t— b;) —Pi— ei), \ (26) 


k= ei, 5; , 
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equation (25) can be rewritten in the much more compact form: 


N 
UW) = it wy. (27) 


Physically, the term 


represents the rate of storage in the large arteries or “arterial uptake.” 
The terms 2;, of course, represent rates at which fluid is being drained into 
the smaller arteries. Notice that if the v; are linear functions of P, then as 
all k; = 0, equation (27) reduces to (5) with 


1 1 
ae 


Steady-state conditions and stroke volume. Before discussing explicit forms 
of v;, we shall consider the conditions which must be met in order for the 
system to be in a steady state. These conditions will be of use for the de- 
termination of the stroke volume and possibly other circulatory pa- 
rameters. 

Let T, be the duration of systole, T the duration of the total cycle, and 
Q(t) the cumulative amount of fluid injected into the tree from time ¢t = 0 
to t = t. Then Q(T,) will be the stroke volume and we have 


Q(T.) = fC vd= f(t v) dt. (28) 


In a steady state, P is periodic with period T and thus 


feo ei —P (Ub. | db 0" 


hence 


ut AP; Ts AP; Lis AP; 
if zeta St aE, at fi Es, a ae 


Equation (29) is obvious from physical considerations; it simply states 
that the net storage in the arteries at the end of one cycle is zero. In other 
words, the amount of fluid injected into the system during one cycle must 
flow out by the end of that cycle, and thus 


Q(T.) = feo) a+ fo (S0,) di. (30) 
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Let K be the ratio of systolic to diastolic drainage, Wey, 


Ts 
ip D v;) dt 
0 


K= , (31) 


Then (30) and (31) yield 
3 
T) =(1 D0;) dt. 
O(T,) = ( on a v,) dl (32) 


Equations (28), (30), and then (29) yield the relations 


iP (S,) a= f(s ee )a= =f (Cla di (33) 


and thus from (32) the stroke volume takes the form 


or) =a+K) f' (es) a. (34) 


Estimation of systolic arterial storage. From equation (34) we see that 
the stroke volume is given by the product of a function of the capillary 
drainage (1 + K) and the arterial storage during systole. It is profitable to 
inquire at this point whether the arterial storage is markedly affected by 
the type of arterial branching assumed in the model. More specifically, 
it is of interest to determine whether or not, at least for the purpose of 
arterial storage calculations, all of the major branches may be lumped 
into an equivalent single branch. Aside from the obvious computational 
advantages gained, this question is important in those cases where no 
distensibilities or arrival times have been tabulated, because then the 
number of undetermined parameters would be considerably reduced. This 
is particularly true in the case of humans, where because of the great in- 
dividual variations in distensibilities, tabulated data of average distensi- 
bilities frequently yields erroneous results when applied to a single person 
(Remington et al., 1949). The single tube model was actually applied to 
man by Warner et al. (Joc. cit.) and it would seem desirable to get some 
idea of the order of magnitude of the error involved in making the single 
tube assumption. 

For this purpose we shall assume the systolic pressure curve is parabolic 
and the diastolic pressure curve decays exponentially. These expressions 
are suggested by the rough considerations leading up to the equations (7) 
and (9). However they may be assumed independently of those considera- 
tions for they do capture the gross characteristics of most pressure con- 
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tours. The fact that the assumption of a non-distorted pressure wave in 
many cases does give good results suggests that the finer details of the 
pressure contour are unimportant. Instead of using (7) per se we shall 
write it in terms of the so-called pulse pressure (i.e., the maximum minus 
the minimum pressure) 6P and the time T,, at which the maximum or 
“systolic” pressure occurs. Thus we have 


Dior 6P 
= Py a f2 OST Gs 
Ey a La iis a 
and 
P(t) =c+ (P,— c) eWETs) In SSI 


Notice that the parameters in the above equations are not all independ- 
ent. That is, the P, in the latter [(9)] can be obtained from (35) by setting 
i = T,. If we use this relation together with the condition Tn < T;, we 
may arrive at the following expression for T,, 


Le ee zip (36) 


To compute the arterial storage we will use (9), (35), and the steady-state 
condition P(T) = P(0). We also have at our disposal the fact that the 
orders of magnitude of ye; and yd; are sufficiently small (e.g., the order of 
7 is around 1 sec. and the orders of the e; and 0; are of the order of 0.1 
sec.) to make the approximations e’: = 1+ ye; and e®% =1-+ yd. 
Carrying out this computation, the following expression for the arterial 
storage during systole is obtained 
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where, alternatively, the value of P, can be obtained from (35), or the 
value of T,, is given by (36). 

If the volume of the equivalent. single branch is equal to the combined 
volumes of the multi-branch model, then letting 1/E be the distensibility 
of the single branch we would have 
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The arrival time ¢ at the end of the equivalent tube (i.e., the time re- 
quired for the pulse wave to travel through the equivalent tube) is not as 
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easily determined as E. We shall assume that ¢ is equal to a weighted su 

of the arrival times to the end of each terminal branch in the multi-branch 
model. The weighting factor in this sum for the ith terminal arrival time 
will be given by h;, the fraction of the total drainage that drains out of 
the zth branch. This is an assumption similar to that made by Remington 
and Hamilton (1947) when all of the capillary beds were lumped together, 
and values for this quantity are tabulated in their paper. The appropriate 
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Ficure 2. Schematic diagram of model similar to that proposed by Remington and 
Hamilton (Joc. cit.). 


LABILE, 1+ 
ARCH THORAX ABDOMEN HEAD VISCERA LEGs 
(1/E) {cc./[@mm. Hg) ]m?} 16 ails) 04 .20 05 .09 
b (sec.) .000 .010 041 .010 O41 .063 
e (sec.) .010 .041 .063 .046 .067 .105 
h 0 0 0 43 .43 14 


* Data taken from Remington and Hamilton (Joc. cit.). 


expression for the storage in the equivalent tube may be obtained from 
(37) by letting V = 1, b; = 0, e; = yg, and 1/E; = 1/E. Then 
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In order to compare (37) and (39) let us consider the multi-branched 
model proposed by Remington and Hamilton. Briefly this model consists 
of the six branches indicated in the schematic diagram of Figure 2. Per- 
tinent values of the parameters calculated from the data of Remington 
and Hamilton (1947) are tabulated in Table I. The 1/Z; are average values 
which have been divided by the square meter surface area of the dogs 
from which they have been taken, and thus cardiac output calculations 
involving the E; are in terms of stroke index rather than stroke volume. 
They are generally applicable within a pressure range of 40-140 mm. Hg. 
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The tabulated arrival times were taken from measurements made at a 
diastolic level of 80 mm. Hg. The arrival time at the end of the equivalent 
tube was computed by assuming 


g= zh; 


where /; represents the relative drainage of the ith capillary bed, during 
one cycle, Le., 
T Vv; 
h; => ae ai dt z 


(See Table I for the assumed values of #; which were taken from the 
Remington-Hamilton papers.) 

Assuming a diastolic pressure of 80 mm. Hg, and using (37), (39), and 
the data tabulated above we have the following expression for the per- 
centage error involved in lumping the six branches into one for the pur- 
poses of arterial storage calculations 
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Since T,, and T,,, are of the order of .1 or .2 sec., it seems apparent that 
with plausible values of P, and 6P the first two terms in the denominator 
will usually dominate any other term in the fraction and yield a small 
percentage error. For example, for a typical case where P, = 110 mm. 
Hg, 6P = 40 mm. Hg, and 7, = .2 sec, the value of T,, calculated by 
(36) is about 0.133 sec. If these values are substituted into (40), the per- 
centage error is of the order of 2%. It seems likely that most other 
plausible examples will give an error of similar magnitude and hence that 
the single tube model represents a good approximation. In any particular 
case, of course, equations (37) and (39) can be combined with the data of 
Remington and Hamilton (Joc. cit.) to yield an estimate of the error. 

Since, as we have seen in many cases, the branching of the tubes has 
only a small effect on the arterial storage calculations, it is desirable to 
determine how sensitive the single tube expression for the storage is to 
errors in the estimation of g. For example, knowledge of this sensitivity 
would have direct bearing on the validity of the pulse contour method of 
H. C. Bazett et al. (1935), which uses the finite pulse wave velocity to 
estimate arterial distensibilities but effectively sets ¢ = 0 in the arterial 
storage and stroke volume computations, 
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If W is the error in estimating y, and Q is the consequent error in the 
calculation of arterial storage, then Q is given by 
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The weighted transmission times tabulated by Remington and Hamilton 
(Joc. cit.) range from .084 to .029 msec. over a pressure range of from 60 
to 180 mm. Hg. If, for example, we took a fixed average value, .056 sec., 
for all time delays of the equivalent tube, then when the “‘true”’ value of 
¢ was equal to .029 sec., and in the case (P,; — Po)/5P = 2 and T, = .2 
sec., we would have an error of about +7%. If, in this example, the true 
value was equal to .084, our error would be only —3%. If we completely 
neglected the velocity of the wave, while ¢ = .063 sec. (as calculated from 
the data in Table I), the error would be about 18%. 

We may conclude that in the case of a non-distorted wave the arterial 
storage can frequently be calculated from the single equivalent tube and 
that this calculation is insensitive to errors in the estimation of the equiva- 
lent transmission time. Further, a crude estimate of the arterial storage 
can be obtained by neglecting both the branching and the finite wave 
velocity. In any particular case, an estimate of the error can be obtained 
by the use of equations (37), (39), and (41) and the tabulated data of Rem- 
ington and Hamilton. Notice that these conclusions apply to the steady 
state, but are independent of any assumptions regarding arteriolar and 
capillary drainage. 

The relation between arterial pressure and arteriolar drainage. In 1933 
Whittaker and Winton perfused the hind limb of the dog with constant 
pressures and found that within physiological ranges these pressures were 
a linear function of the ensuing flow. If the linear portion of this curve was 
extrapolated backward, it intercepted the pressure axis at about 20 mm. 
Hg. Citing this experiment, or relying on Poiseuille’s law, several authors 
(Bazett et al., 1935; Remington and Hamilton, Joc. cit.; Warner et al., 
loc. cit.) have assumed the result equally valid for pulsating pressures and 
hence set 


W=E [P (t— $:) — 4] (42) 


where R; is a constant and P(t — ¢,) is the pressure at the head of the 7th 
arteriolar bed, ¢; the arrival time to that bed, and @ is a constant of the 
order of 20 mm. Hg. 
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Let 1/R = = 1/R;and let D be the total amount of fluid draining into 
the arterioles during one cycle. Then in a steady state P will be periodic, 
and hence 


a il : 1 u 
D= zg f PU) — a a zi PO-4 
If %; is the fraction of fluid draining into the 7th bed, 
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and hence the expression for the velocity of cardiac ejection can be 
written as 
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The unknown XK in the expressions for Q(T.) is then given by 


shi f' (Pt #) — 6) dt 
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(45) 
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The constant R can be determined from equations (33), (42), and (43), 
Le), 


fi Bh) — 1 ai 
R= , : (46) 
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If we are given the pressure contour, then equation (44) together with 
(46) would be sufficient to determine the cardiac ejection curve, and equa- 
tions (45) and (34) determine the stroke index. In the following para- 
graphs we shall develop expressions for these quantities which are easier 
to deal with, and perhaps even more general in their range of applicability. 

Performing power series expansions of AP; and P(t — ¢;) around #, 
equation (44) can be rewritten as 


(je yu, eared (47) 
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where 
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Since all of the arrival times are of the order of 0.1 sec. or less, this series 
probably converges very rapidly. For example, using the values in Table 
I, and taking R of the plausible order of 1 (mm. Hg)(sec.)(cc.—), the ratio 
of a3 to az is of the order of .0005. Hence if the higher order derivatives are 
not too large, it would be reasonable to take only the first two terms of 
the summation in (47) and write 


d?P = Moco ty) 
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At the beginning of systole there is a sudden rise in pressure which implies 
that the higher order derivatives are probably very large, and thus during 
early systole (48) would not be applicable. However, in diastole, shortly 
after the dicrotic notch the pressure changes relatively slowly, frequently 
decaying with a time constant of the order of 1 second. During this period, 
all derivatives would be somewhat of the same order of magnitude and 
equation (48) should closely approximate the situation. Since U = 0 
during diastole, the solution to (48) is of the form 
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where 7, m, My, and WM are constants. 

Using plausible orders of magnitude, it is apparent that a; and 1/R 
generally will be positive and a2 negative. Hence the roots m and m, of 
the characteristic equation of (48) will be real, one of them, mz, must be 
positive, the other, m, must be negative. Therefore, in order to avoid the 
possibility of infinite pressures, the coefficient M, must be zero. 

Suppose a pressure contour is measured and an exponential decay of 
the form (49) with M, = 0 has been fitted to its diastolic portion. Then, 
substituting (49) into (48) and solving the ensuing expression for R, we 
arrive at 


pees! te (50) 


As would be expected, R will always be positive because m is negative. 
The value of m is obtained from the pressure contour, and all of the other 
quantities on the right-hand side of (50) are obtainable from the pub- 
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lished data of Remington and Hamilton (Joc. cit.). Hence R is completely 
determined. The fitting of (49) also enables us to arrive at a value for 8, 
and since in a steady state Q(T.) is equal to the total amount of fluid 
draining out of the arterial tree, we have 


or) =2>>% afte Peo) 0) d= gf PoO-sas (51) 


as an expression for the stroke volume. 

Note that in this determination @ is obtained directly from the pressure 
contour, but in the previous formulae 6 was assumed to be known. Based 
on the experiments of S. R. F. Whittaker and F. R. Winton (1933) several 
authors (Bazett et al., 1935; Remington and Hamilton, loc. cit.; Warner et 
al., loc. cit.) have assumed 6 = 20 mm. Hg in all cases. However, the per- 
fusion experiments of J. R. Pappenheimer and J. P. Maes (1942) reveal 
that extrapolation of the linear portion of the pressure flow curve yields 
an intercept on the pressure axis which increases with vasoconstriction. 
Pressure intercepts higher than 100 mm. Hg have been found. They sug- 
gest that the Whittaker-Winton experiments were performed on strongly 
vasodilated preparations. (In this connection see also Folkow, 19520.) 
Hence, in general, we may expect @ to be higher than 20 mm. Hg. This is 
confirmed by the measurements of Apéria (Joc. cit.) who finds that @ as 
measured by fitting an equation of the form (49) (with M, = 0) usually 
ranges between 60 and 80 mm. Hg. The present approach is thus advan- 
tageous in that it does not specify a priori the numerical value of 6 and 
that it allows both 6 and R to vary under different physiological condi- 
tions. Further, this approach is not restricted to the steady state. That is, 
for the non-steady state, assuming R remains relatively constant within 
a single cycle, a combination of (50) and (44) yields an expression free 
from undetermined parameters for the rate of cardiac ejection, and its 
integral would then give the cardiac ejection curve. 

Just as in the case of arterial storage, it is of interest to determine the 
discrepancy in calculating R between the multi-branched system and a 
single-tube equivalent system with distensibility 1/E = D 1/H;, and 
transmission time g = Lhid;. If Rm is the R calculated for the multi- 
branched system, and R, is the R calculated for the single equivalent 
branch, then from (50), 
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Using the values in Table I and letting m = —2 sec., for example, the 
error, as calculated by (52), is less than 1%. A value of m = —2 sec. isa 
rather fast decay, but for slower decays the error given in (52) would be 
even smaller. Thus in so far as stroke volume calculations are concerned, 
there seems to be little advantage in preferring the multi-branched system 
over the equivalent tube. For a single tube (50) reduces to 
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The values of 4; are assumed to be known. Further, the specification of 
the ¢; in terms of the transmission times through the various branches 
is somewhat arbitrary; (cf. Remington and Hamilton, Joc. cit.. Warner 
et al., loc. cit.). For these reasons it is of interest to determine how sensi- 
tive the expression for R is to errors in the estimation of h;, ¢;, and g. Let 
mg — 1 be denoted by w, and again let y be the error involved in estimat- 
ing ¢, then the error in calculating R will be given by the expression 
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If we consider, for example, an error of — 100% (in other words letting the 


estimated g = 0 and hence y = —1), expression (54) reduces to 
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Note that because m < 0, w is always negative. Direct calculation shows 
that the absolute error computed from (55) will increase monotonically as 
w decreases from —1 to — ~. Again using the values in Table I and as- 
suming a rather fast decay by setting m = —2 sec.—', the error given by 
(55) will be only about 6%. For slower decays w will be closer to —1, and 
consequently the error will be even less than 6%. 

Contingent upon the validity of equation (48), it would therefore ap- 
pear that in so far as stroke volume determinations are concerned, the in- 
troduction of a branching system with finite pulse wave velocities in many 
cases differs from a lumped system by at most 10%, provided the distor- 
tions in the pressure wave are neglected. For stroke volume computations, 
the important consequence of a finite wave velocity then would not be in 
the fact that there is a time delay between systole and the appearance of 
the pressure pulse (at a particular point down stream) per se, but in the 
possibility of such a situation leading to wave distortions through local 
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reflections (Porjé, 1946; Karreman, 1952, 1953, 1954). However, we must 
also face the fact that a linear pressure flow relation in the arterioles is 
not completely justified by an allusion to Poiseuille’s law, or the citation 
of the experiments of Whittaker and Winton, as has frequently been as- 
sumed by many authors. 

The distensibility of the arterioles together with the complications of 
an anomalous viscosity and the existence of critical closing pressures 
(Burton, 1951) suggest that the steady flow should be a non-linear func- 
tion of the pressure. Green et al. (1944) and Folkow (1952a) have fit their 
data by relating the flow to a fractional power of the pressure. The frac- 
tional power according to Folkow can have much more variation than 
previously thought. This is not too distressing since it still may be argued 
that the pressure flow relation may be reasonably linear over the pressure 
range occurring in one cycle, and since we have not rigidly fixed the values 
of either R or 0, the possible variation in the parameters of the approxi- 
mate relation v = (P — 6)/R is permitted. 

However, the fact remains that Poiseuille’s law is only valid under 
steady pressures, and the experiments cited above were carried out under 
steady pressures; yet we are dealing with pulsatile pressures. The forces 
involved in propelling the blood through the capillary beds consist of a 
driving force (determined by the pressure gradients) and a viscous re- 
tarding force (determined by the velocities of flow). Under steady pres- 
sures, these two forces eventually balance each other and a ‘“‘Poiseuille” 
type flow ensues. Under pulsatile pressures they do not necessarily balance 
and hence the inertia of the blood may become important. Green et al. 
(Joc. cit.) have studied pressure flow relations in the perfused hind limbs 
of dogs. He found, over a long period (1.e., several cycles), little difference 
between the cumulative amount of fluid draining into the capillaries under 
pulsatile pressure and the amount draining in under steady pressure. From 
this result it might be suspected that under normal pulsatile conditions the 
viscous forces are completely dominant and hence it is justifiable to apply 
the steady-state solution. However, we shall soon see that this is not the 
case. 

In order to more clearly discuss the various possibilities, let us consider 
a simple case by assuming the system can be treated in a “lumped” fash- 
ion, and define P; as the “pressure in the arterioles.” Then if P is the 
pressure at the terminal end of the artery, we have 
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where A and B are assumed to be approximately constant over any par- 
ticular cycle, and v denotes the fluid velocity into the bed. It should be 
pointed out that in this crude treatment linear velocity (cm./sec.) is 
assumed proportional to volumetric velocity (cc./sec.). 


Under steady pressures, and after sufficient time has elapsed, (56) re- 
duces to the familiar results 


v= (P-P). (57) 


Here, the pressure intercept of the extrapolated linear pressure flow curve 
is given by P;. There is no reason to expect P; to be constant (e.g., equal 
to 20 mm. Hg) under all conditions, in fact the experiments of Pappen- 
heimer and Maes (Joc. cit.) require that P; rises during vasoconstriction. 

Referring to the results of Green et al., consider the differences in cumu- 
lative flow if the bed is perfused with steady pressures or with pulsatile 
pressures over a period of several cycles. In the case of steady pressures, 
by integrating equation (57) from 0 to wT + ¢ where x is an integer and 
Oe << 7, we have 
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where P and P are the steady perfusion pressures which equal the means 
over one cycle of corresponding pulsatile pressures. The same integration 
performed on (56) yields 
ar +t d v (ite nT +f 
es = f= dt 59 
i aaAf | @-P)d Bf (59) 
or, since in a steady state, both P and v are periodic, 
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If n is sufficiently large, the differences between (58) and (60) become 
negligible. Thus we see that the result of the aforementioned experiments 
of Green ef al. can be anticipated without any assumption as to the rela- 
tive magnitude of the viscous forces. 

We are concerned with pulsating pressures in a very complicated net- 
work of fine tubules. However, according to B. W. Zweifach (1949), by the 
time the blood reaches the narrow arterioles pulsatile flow has practically 
ceased. This suggests that the response (in terms of fluid flow) of the ar- 
teriolar-capillary bed to change in pressure is sufficiently slow so that the 
pressure pulsations in the arteries induced by the intermittent beat of the 
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heart are “imperceptible” to the bed. In other words, rather than respond- 
ing to the pressure at any given moment, the drainage responds to some 
function of the mean pressure Pn. (Credence may be lent to this interpre- 
tation by the experiments of Williams and Schroeder, 1948.) These 
authors measured the pressure decay rate distal to a suddenly imposed 
arterial constriction. They found the instantaneous pressure at the mo- 
ment of occlusion has a less significant effect on the pressure decay curve 
than the mean pressure taken over an interval just prior to the occlusion. 

If we do assume that the drainage is some function of the mean pressure 
Pm, more specifically if 

0;= 0;(Pm) » (61) 

where v;(P.) is some unspecified function, then it will be seen in what is 
to follow that no specification of ;(P,) or Pm is necessary, and hence we 
need not concern ourselves with any of the complications mentioned 
above. Further, knowledge of the relative drainages of the various beds 
(the /;) will not be necessary. 

That is, we may combine (61) and (28) to form 


je Udt = fa (Sat D (Pm) Ts (62) 
but, from (30), 


fua- v(P,)T. (63) 
Solving for v(Pm) in (62) and (63), we find 
Pn) =p f (Sepa, - (64) 
and our expression for the ejection velocity becomes 
U = Lex Ca Sf ‘aP.at). (65) 


Finally, in this case the value of K, the unknown in the expression (34) 
for Q(T), can be obtained from (31) in the form 


Ts 
i iit (66) 


In this formulation we have ignored the events taking place between the 
terminations of the major branches and the terminations of the arterioles 
where there is no pulsatile flow. It is primarily this problem that motivates 
the following section which attempts to remove at least some of the re- 
strictions imposed by the assumption of a non-distorted wave, and at the 
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same time take into account the distensibility of the smaller arteries and 
arterioles. 

Damped and distorted pressure waves. As the pressure wave travels 
through the network, its form is distorted by local reflections, changes in 
distensibility, and pulse velocities. Further, the smaller arteries offer con- 
siderable frictional resistance and hence there must be a viscous damping 
of the wave until finally, according to Zweifach (1949), the pulsations are 
no longer present at the ends of the arterioles and beginning of the capil- 
laries. It is at this point that the outflow (in a steady state) is independ- 
ent of time, and it is here that we shall assume our network of tubes ends. 
Thus the terminal outflow velocities v; are considered to be constant. 
Since the branching is highly complex and unspecified, we have no way, 
even in principle, of calculating the pressure in each segment. However, 
we can measure the pressure contour at the origin of the tubes and, pre- 
sumably, the pressure contours in the various segments bear some relation 
to the original contour. 

Let us fix our attention for the moment on all the segments which have 
the same arrival time r. These segments all will be very roughly the same 
distance from the origin and hence it is within reason to assume that they 
have very similar characteristics. In particular we shall consider the case 
where the £;; are functions of 7 alone, and if the pressure contours of all 
segments with arrival time 7 are summed, the distortions tend to cancel 
out. This last condition introduces an element of randomness and sums 
up our lack of knowledge of the intimate details by stating that, if we are 
likely to have a certain distortion in a particular segment, we are just as 
likely to have an equal and opposite distortion in another segment which 
has the same arrival time. In other words, we are assuming that, if the 
pressure contours in all segments with the same arrival time are added, a 
damped version of the original contour is obtained. However, we may re- 
lax the restrictions a bit further by allowing a time independent, but 
otherwise arbitrary, additive weighing factor (which may vary from seg- 
ment to segment) to be added to the pressure in each segment before the 
‘above summation is carried out. If we let # be the pressure at the terminal 
ends where no pulsations occur, and if P(¢) = p(t) + ? denotes the origi- 


7) 


nal pressure contour, then by introducing the symbol Ss to denote sum- 
mation over all segments with arrival time 7, we may write 


(7) 
D> Paty = elt) op (t—7) +h (67) 
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where g(r) is simply a damping factor, and 4; is the additive weight for 
the ij segment. Letting f(r) = [g(7)/E:;(7)] we obtain 
DEEP ed 
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Since we are assuming the validity of (68), it will be worth while to indi- 
cate a slightly different set of conditions which gives rise to the same ex- 
pression. That is, suppose the Z;; are not necessarily functions of 7 alone, 
and let P‘7 represent the pressure in any segment with arrival time r. 
Then we may define the function z;;(¢) by the relation 


PY) = guP @— 7) +24 (69) 
where gi; is the damping factor of the ij segment. We will consider the 


case where the z;;(¢) are nearly constant, or more generally where 
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(x) 
and identifying }° g:;/E,; in (70) with f(z) in (68) reveals the equivalence 
of the two relations. 

At this point let us reintroduce equation (15) and recall that it is es- 
sentially an equation of continuity which must hold no matter how vio- 
lently the pressure wave may be distorted in the various segments. It re- 
quires a summation of the pressure derivative (weighted by the proper 
factor) over every volume element. The order in which this summation 
is carried out is immateral. In particular, we choose to first sum over all 
volume elements which have the same arrival time and then over all ar- 
rival times. Carrying out this procedure with the aid of (68) we have 


U= Sp OZPU-1) $09 (71) 


(x) 
where u(r) = S* 0. 

Expression (71) may be converted into an integral by the following 
considerations: If Ax;; is the length of the zj segment and Ar;; is the time 
required for the pulse wave to pass through this segment with velocity 
Cis, then 

AX ig = CiyjAT i - 
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We may define the factors p;; and Ar by the relation 
Dae = pijAt 


where the value of Ar is independent of i, but as all Ax;;— 0, Ar > 0 
(e.g., Av could be defined as the average of all Ar;,). Combining the last 
two expressions yields 

Ax j= Crp QeINGe « (72) 


Letting 1/;; be the distensibility per unit length of the uy segment, we 
also have (by definition) 


Ei; = Ax sj 3 

At this point we must consider separately the situations leading to (68) 
and to (70). In the former case, if in passing to the limit we specify that 
all Ax;; are equal and hence independent of 7, it follows from (72) that the 
product c;;o;; must also be independent of 7. Then letting 


g(r) 
F(z) =e. 
it can be seen that 
f(r) =F (rn) Ar. (73) 


In the situation leading to (70) we can arrive at the same form as (73) by 
simply identifying F(r) with 
(r) 
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Inserting (73) into (71) and passing to the limit (1.e., letting Ax;; — 0) 
we have 
U() -U= [PO GPU d (74) 


where 7 is the maximum 7 and U = >) v(r). It is apparent from physical 


considerations that in a steady state U must equal the mean value of U. 
This can be readily verified by integrating (74) from 0 to T, the duration 
of one heart cycle. Equation (74) is a Fredholm integral equation of the 
second kind with the kernel d/dé[P(¢ — 7)]. If P(@) and U(¢) are meas- 
ured, (74) can be solved by numerical methods for the unknown. function 
F(r). Knowing the shape of F(r) and how it varied with different condi- 
tions would possibly be of aid in the further development of the theory. 
This function is determined not only by the distensibility of the arterial 
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tree but also by the capacity of the tree to cause distortions through re- 
flections and viscous damping. 

By definition, F(r) > 0 throughout the interval 0 < r < 7 so that we 
may apply the mean value theorem to (74) and arrive at 


UW) —0= FP =) [UF (1) ar (75) 


where »; is a mean value of 7, 1.e., 
0 < Vt < nH - (7 6) 


Note that », is a function of ¢. That is, we must perform the integration 
over 7 for each value of ¢, and hence we must anticipate a dependence of 
the extracted mean value on ¢. Let us define the constant po by the 
relation 


7 
= / F at; 
Ko i (r) dr 
then (74) may be rewritten as 


UW) —0=m SP»). (77) 
Integrating (77) from 0 to ¢ where 0 < ¢ < T we have 
t 
J (U-D) dt= wo (Pi—») —P(—») 1 . (78) 


Equation (78) states that the integral of the excess of the instantaneous 
fluid velocity over its mean is a linear function of the pressure at a previ- 
ous time. Thus we would anticipate this integral curve to have the same 
form as the pressure curve but lag behind it by an amount 1. 

Figure 3 is a reproduction of some published curves of R. E. Shipley, 
D. E. Gregg, and E. F. Schroeder (1943) taken from the femoral artery 
in the dog. It is seen that despite the apparent complexity of the velocity 
curve its integral follows the pressure curve. Further, there is a time de- 
pendence of the lag between the two curves, for example, the peak of the 
flow curve lies a measurable amount behind the corresponding pressure 
peak, whereas near the end of the cycle the time lag is practically zero. 
Thus, at least in this instance, our model appears to adequately describe 
the situation. 

We have not calculated the dependence of »,, the mean value of r ex- 
tracted from the integral in (74) on ¢; however, the variation in »; is greatly 
restricted by inequality (76). 

The propensities of the solution of (74) can be more easily studied by 


PRESSURE AND FLOW IN THE CIRCULATION 197 


the following means. Denoting dP/dt by P we may expand the function 
P(t — r) in a power series in terms of 7 around the point (r = x) where x 
is an unspecified constant. If this is done, (74) can be rewritten as 


U-U= > fro Se or (St) @ 


or, since 


is independent of r, 
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Ficure 3. Femoral artery pressure flow patterns taken from paper of Shipley e¢ al. (loc. 
cit.) with the authors’ and publisher’s kind permission. 


Defining yu, by the relation 


o/] 
= eh) Ra & 5 GS) 
pin = f° (1 — 6) °F (1) der 

noting that 

d . dE 

BP-n =-Fton), 
and remembering the definition of P, we may rewrite the above expression 
as 


oe) 


De ni oP make (80) 


n=0 
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If we choose x such that py = 0, then, by (79), «x will equal the mean value 
of r weighted by F(r) over the range 0 < 7 < 7. It follows that «<9; 
further 


0 
me d 
| oe °F (7) Kid 
Ho} q 
F d 
ue (r) dr 


and, therefore, if 7 is small, the series (80) must converge very rapidly. 
As seen in Table I, the maximum transmission time between the root of 
the aorta and the terminations of the major branches is of the order of 
1 sec. The distance from this termination point of non-pulsatile flow is 
very short compared to its distance to the aortic root and hence unless 
the velocity of pulse-wave transmission suddenly becomes very small in 
the finer branches, we would anticipate an 7 close to .1 sec. (if the origin 
of the tubes is taken to be the aortic root). Further, the closer the origin 
of the tubes to the capillary beds, the smaller the 7. Bazett and Dreyer 
(1922) have found that velocities in the larger vessels are slower than 
those in the smaller peripheral vessels, but in the more peripheral parts 
the velocities were more variable depending on conditions of local vaso- 
dilitation. Thus it is rather unlikely that the above objection (a sudden 
showing of the pulse wave in smaller vessels) to our estimation of the 
order of magnitudes of 7 holds. 

The ratio of the fourth term to the first term in the series (80) is then 
of the order 73/6 ~ .0002 and possibly even less. Hence it would appear 
reasonable to take the first three terms of the series as a first approxima- 
tion and write 


nn (81) 


OPA =tk) 4 Me a®P (t— x) 
dt? 2 dé 


= d 
U-U= oa PU—«) +m (82) 
Since we have chosen « in such a manner as to make p, = 0, equation 
(82) has the same form as (75), except that the ambiguity in the time lag 
has been removed (.e., « is independent of #) by adding a small correction 
term of the form 


If the correction term is ignored, and equation (82) is then fit to the data 
in Figure 3, the order of magnitude of « is roughly .02T and yo ~ .02 cc./ 
mm. Hg. These values are for that part of the vascular tree lying between 
the point on the femoral artery where the measurements were made and 
the point where non-pulsatile flow occurs. 
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If we choose x = 0, then of course 4, + 0, but the above considerations 
with respect to the rate of convergence of the series still hold. We are, of 
course, at liberty to take the origin of our tubular system at any point. In 
particular by taking the origin of the tubes closer and closer to the capil- 
lary beds, » becomes smaller, and the series converges faster. This sug- 
gests that rather than using a Poiseuille type flow as an approximation 
for the boundary conditions in arteriolar drainage, the simple expression 


U (i) -0=mo (83) 


should be employed, where now U, U, and P represent the instantaneous 
flow, mean flow, and pressure respectively at the head of the arteriolar 
bed. It is possible that the order of magnitude of uo(dP)/(dt) may be very 
small when compared with U(t) — U. In such an instance (83) would re- 
duce to a previously considered case [i.e., eq. (61)]. 

The heart parameters and their relation to the arterial pressure. Return- 
ing for a moment to conditions at the energy source, let us recall equation 
(1), where it was assumed that the heart’s radius 7 diminished at a con- 
stant rate \ during systole. Because the cardiac wall is muscular, it is to 
be expected that when the initial conditions (e.g., the diastolic volume) 
are altered, the value that \ assumes during the contraction will change. 
As indicated in the introduction, it would be desirable to have a rough 
idea of the dependence of ) not only on the pressure in the arterial reser- 
voir but also on the initial conditions of the muscular wall of the heart. 

With this end in mind, let us apply the velocity load equation of Hill - 
(1938, 1949). That is, if Z is the average tension (load) on the fibers in the 
cardiac wall, and if Lo is the tension that would be developed by the fibers 
if they contracted isometrically, then Hill’s equation states that the con- 
stant velocity \ of an isotonic contraction will be given by 


(a+ZL) r= 6b (L)—-L) (84) 


where a is the constant heat of shortening, and 0 is also a constant. 
This equation has been tested in various species and found valid both 
for whole muscles or isolated muscle fibers (Wilkie, 1954). The experi- 
ments of G. Lundin (1944) indicate that for cardiac muscle Lo varies 
linearly with resting length over a considerable range, and hence we shall 


let 
Ly) = Bro—o (85) 
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where 6 and o are constants. If the heart contracts as a spheriod, and the 
interventricular pressure is denoted by P., then L will be given by 


L=aPy,r (86) 


where a is a constant (equal to 4 for a perfect thin-walled sphere). 
Eliminating Z and Ly from (84), (85), and (86), and solving for 

yields 

ne b(Bro— o — aPyr) 


7 
a+ aPyr ea) 


r 


In particular, at t = 0, P, is equal to Po, the diastolic aortic pressure 
(more precisely, the valves open the moment P, just exceeds Po) andr = 
ro. Substituting these values into (87) we arrive at an expression relating 
solely to the initial conditions; in other words, we have the velocity of 
contraction as the following function of Po and 70: 


b (Ce) / Aime gee 
ro a a 


fo Pi rg 
a 


r (88) 


It is most likely that (84) and hence (87) do not hold at exactly ¢ = 0; 
however, we may assume that they do hold at a slightly later time when 
the product P,r is still sufficiently close to Povo. A combination of (4) and 
(88) gives an expression for the cardiac ejection curve in terms of the 
heart’s parameters, the diastolic arterial pressure, and the diastolic radius 
of the heart. These relations can be checked experimentally, and, if suc- 
cessful, afford a means of estimating heart parameters which have mean- 
ingful interpretations in terms of the actual muscular elements of the 
heart. 

As stated above, in order for the valves to open and for any blood to 
flow, P, at t = 0 must exceed the arterial pressure. It is of interest in 
passing to note that this condition implies an inequality among the heart 
parameters which may be of value in future applications. That is, if any 
shortening is to occur, equation (84) requires that the maximum possible 
isometric tension Ly be greater than the tension Z actually existing in the 
fibers. Thus using (85) and (86) we have 


Ilo=Bro—o=L= aP,7ro 


or 
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For the valves to open 
dhol ane 


F-(2)>p, (89) 


which relates the parameters of the heart to the pressure against which it 
may pump. 

Thus, for a given ro, there exists a critical pressure P;*, such that if 
P,) > Ps, no blood will flow from the heart to the arteries. By plotting 
corresponding values of 1/r) and P;, a rough estimate of o/a and B/a may 
be obtained. 

In previous sections we have seen that simple considerations regarding 
the heart, when combined with the air chamber theory, give rise to a 
pressure contour which roughly approximates the qualitative shape of 
those frequently found in experimental measurements. However, the sig- 
nificance of these equations is greatly limited until their parameters have 
been evaluated. We shall utilize the results of previous discussions in the 
evaluation of these parameters by setting the stroke volume predicted by 
the air chamber equations equal to the stroke volume predicted by the 
pulse contour methods. In other words, we are assuming that the pressure 
contour obtained from the air chamber theory will have some rough 
quantitative significance if it corresponds to plausible flow values. 

The two situations we are most likely to be confronted with are: 
(1) where the parameters of the heart are given and the pressure contour 
sought; or, conversely, (2) where the pressure contour is given and the 
heart parameters sought. 

Considering first the former case, suppose the diastolic heart radius 70, 
the velocity constant \, and the diastolic pressure Py are known. Note 
that if the cardiac muscle parameters, a, b, a, 8, and o are known, \ may 
be calculated from equation (88). In addition, the magnitudes of the total 
cycle duration T and systolic duration 7, are assumed known. If T,, is not 
known, it probably can be estimated from one of the many empirical for- 
mulae which have been proposed to relate T and 7, (Katz, 1940). For a 
perfect sphere, A may be set equal to 47. Then the only remaining parame- 
ters of equations (7) and (9) are v, ¢, Ps, and c. If the value of ¢ is as- 
sumed known, we will be able to evaluate the other three unknowns by 
making use of the pulse contour method. Despite theoretical objections 
raised in previous sections, these methods have in many instances given 
reliable results and hence are worthy of consideration. 

If the theoretical contour is given by equations (7) and (9), then appli- 


Therefore we must have 
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cation of the pulse contour methods would yield a value Qp for stroke vol- 
ume which is a function of the variables y, e, and P as well as the tabu- 
lated distensibilities and arrival times. On the other hand, the stroke vol- 
ume of our spherical heart is given by 


3 (ri (Ha aDy, 


Then one constraint on the unknowns ¥, e, and P, would be 
Or (ys 6 Pa) = 3 [AS— (ro— 07) 4] . (90) 


A second independent constraint would be the steady-state condition 
P= PO vorsirom (9) 

(P,— c) e-¥?-7T) + c=Py. (91) 
Finally, the continuity of P(é) yields the third independent constraint, i.e., 


P,=Pot [eAdro— y (Po— 6) ITs (92) 


—4 [yeArg— 7 (Po— c) + 2eAN rol T: 


Thus the parameters y, e, and P, can be solved for and hence the whole 
contour is determined. (Similar considerations can be applied to the case 
of a cylindrical chamber.) 

The second situation mentioned above assumes that the contour P(E) 
is given, and values of the heart parameters are unknown. In this case, 
P., Po, T;, and T can be read directly from the pressure curve, and hence 
y and ¢ can be determined by fitting (9) to the diastolic pressure decay. 
Then W may be evaluated by use of the stroke volume Q> predicted by 
the pulse contour method. Specifically 


1 (P—c)dt (93) 


Since e = y/W, we have now determined ¢, y, and W. Next a parabola 
is fitted to the systolic portion of P(t). This will be of the form 


P(t) =Pot qit— gol? (94) 


where q and q, are constants determined by the fitting process. If, mere- 
ly for the sake of simplicity, we assume that the parabola is fitted in such 
a way that it passes through the maximum pressure, we may write 


OL Bets 
ss a hi he 
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Comparing (94) to (7) reveals that 
2 
Adri (Pye) = qi = OP (95) 
and 
4 [yeAArp— (Po— c) +2eAN ry] = Wee (96) 
Solving (95) for \ in terms of 7 yields 


6P 
fink q+ ¥ Po c) 


im eA re (97) 
Using (97) in (96) and solving for ro yields the relation 
6P z 
[2+ 0) 
roe=3 PA (98) 
yea a) 


Corresponding equations for a cylindrical chamber would be given by 


As 26P+y(Po- OTn 
Um eN DEP (l= 4T w) 


(99) 


and 


vo V/26P(1— Tn) €A’. (100) 


By recording several pressure contours under different conditions, sev- 
eral sets of values for Po, \, and 7) may be obtained. These sets may be 
then used to evaluate the muscle parameters a, b, a, 8, and o of equation 
(88). 


I would like to thank Professor H. D. Landahl and Dr. Clifford Patlak 
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With reference to several recent papers by the author, it is pointed out that within the 
principle of biotopological mapping a choice of a primordial graph and of a particular trans- 
formation defines a system of abstract biology, similar to systems of abstract geometries. The 
study of such abstract systems is necessary before one can be found which is isomorphic to 
the actual biological world. A brief survey of the structure and properties of the system based 
on the choice of the primordial graph and of the transformation T defined in a previous paper 
(Bull. Math. Biophysics, 16, 317-48, 1954) is made. Two more topological theorems are 
demonstrated, which, interpreted biologically, lead to the conclusion that the higher an organ- 
ism, the more adaptable it is. Finally a criticism of that particular system of abstract biology 
is made, and its inadequacy for the representation of the actual biological phenomena pointed 
out, and a suggestion is made for a possible point set topological approach to biology. 


In a recent paper (Rashevsky, 1954; referred to as I) we pointed out 
that mathematical biology has hitherto emphasized only the quantitative, 
metric aspects of biological phenomena and has entirely neglected the re- 
lational aspects, which are just as important. This neglect of the relational 
aspects of biology results in one-sidedness, which leaves almost untouched 
such problems as the unity of the organism and that of the whole organic 
world. 

It is very essential to know the quantitative aspects of different physi- 
ological or general biological phenomena. But the study of those quantita- 
tive aspects alone does not bring out the important relational aspects 
which connect those single phenomena into a united whole. Such phe- 
nomena as response to stimuli, locomotion, ingestion, digestion, defeca- 
tion, sexual intercourse, etc., are both quantitatively and qualitatively 
very different in different organisms. Yet the basic relations in which those 
individual phenomena stand to each other are fundamentally the same in 
a protozoan, a plant, an insect, a bird, or a man (I, pp. 322-34). 

The natural mathematical tools for the study of such relational aspects 
are topology and group theory. A characteristic property of the biological 
relations is that whereas the corresponding phenomena in different or- 
ganisms may be of very different degrees of complexity, the similarity of 
the relations between different phenomena enables us to establish a many- 
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to-one correspondence of a higher and lower organism with respect to 
these phenomena (I, p. 324). The translation of this well-known biological 
situation into topological language and its combination with the basic 
postulate of every science, namely, that of the existence of uniformities in 
nature, leads us to the establishment of the following general principle, 
which we call the principle of biotopological mapping (I, p. 325): 

If we can in some way represent the different biological relations within 
an organism in terms of a properly chosen topological space or complex, 
then the topological spaces or complexes by which different organisms are 
represented are obtained from one or at most a few “primordial” relative- 
ly simple spaces or complexes by the same geometric transformation, 
which contains one or more parameters, to different values of which cor- 
respond different organisms. 

As we emphasized in I, the above principle is not a hypothesis but 
merely a statement of a known fact coupled with our basic belief in the 
uniformity of nature, without which no science can exist. 

The next problem is to establish within the above principle a set of hy- 
potheses that would form the foundation of a theory. 

In I we suggested a one-dimensional complex, the directed graph, as a 
possible representation of the relations within an organism. The points of 
the graph correspond to the individual biological functions, while the di- 
rected lines show their interrelations. In particular a very hypothetical 
and tentative graph was suggested as a possible primordial one (I, pp. 
325-28). 

A geometric transformation which leads from this primordial graph to 
much more complicated ones has also been suggested, using as a lead again 
a known biological fact, namely, that an increase in complexity of an or- 
ganism is accompanied by specialization of cells and tissues, so that some 
of those cells and tissues lose some of their biological functions and acquire 
some new ones. In terms of the representation of an organism by a graph 
this means that a more complicated organism is represented by a graph 
that is composed of a number of primordial graphs, properly connected, 
and in which some have lost certain points, others have lost other points. 
Considerations of that nature lead to a transformation rule (7) that pro- 
duces from a relatively simple primordial graph a very large but finite 
number of transformed or derived graphs. Inasmuch as any such graph 
describes complicated physiological relations within the organism it rep- 
resents, and inasmuch as those relations can be in principle studied experi- 


mentally, we have basically a theoretical system capable of experimental 
verification. 
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The finite number of possible transformed graphs represents the total 
possible number of different organisms. Expressions for the calculation of 
this number have been derived (Rashevsky, 1955a, b; referred to respec- 
tively as II and III), and it has been estimated at 108 (I, p. 341). This is 
much larger than the total number of known existing species (about 10°). 
Of course, not all possible organisms may have as yet been produced. On 
the other hand, the discrepancy may mean, and possibly does mean, an 
improper choice of the primordial graph and/or of the transformation. 
Right at the outset we thus find an experimentally or observationally 
verifiable consequence of the theory. 

In a subsequent paper (III) some theorems were established in regard 
to the effect of the transformation T on the point base of the primordial 
graph. Making the not unnatural assumption that a biological function 
can be regenerated, if lost, only if the biological functions which lead to it 
remain intact, we arrive at the interpretation of the points of a point base 
as representing such biological functions as are essential to the regenera- 
tion of any other biological function, if the latter is lost through injury. 
Theorem 2 of III leads us then to the conclusion that the ability of regen- 
eration decreases with increasing differentiation and complication of an 
organism, a fact well known in biology. 

The above conclusion depends only on the transformation T, regard- 
less of the choice of the primordial graph. Thus some general biological 
laws may be derived from the study of the transformation only. 

A system of abstract biology. If we choose a different primordial graph, 
or a different transformation, or both, we shall in general be led to differ- 
ent conclusions. Each choice gives us an abstract system of biology, within 
the framework of the general principle of biotopological mapping. Our 
goal is to find a system which is isomorphic, or as nearly isomorphic as 
possible, with the actual biological world. In order to eventually reach this 
goal, it may well be necessary to develop first a general theory of such 
abstract systems of biologies. The system based on the primordial sug- 
gested in I and on the transformation T is just one of the infinite number 
of possible ones, just as Euclidean geometry is only one of an infinite num- 
ber of possible geometries. An abstract study of the systems of geometry 
of necessity preceded their application to the physical universe. A study 
of the abstract systems of biology will possibly have to precede any fruit- 
ful application of topology to biology. 

Let us now more closely examine the structure of the system hitherto 


discussed. . 
The study of the transformation T alone has led us to the general bio- 
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logical law which states that higher organisms have a lower capacity for 
regeneration. We may seek to establish some other general laws, which are 
independent of the choice of the primordial. We now shall prove another 
theorem, the biological interpretation of which leads to another experi- 
mentally established law. Familiarity of the reader with I is assumed. 

Definition. If in a directed graph we have a cycle such that all arrows 
are in the same direction, in other words, a cycle which is a closed directed 
way (III, p. 114), such a cycle is called a uniformly directed cycle. Ex- 
amples of uniformly directed cycles are the cycles ASsD..SaeDA and 
S;CsEM{IDAS;, in Figure 1 of I. In the same figure the cycle F,EM;F, 
is not a uniformly directed cycle. 

Theorem 1. If r specializable points fi, f,...,f, of the primordial 
graph P all belong to a uniformly directed cycle, and if no two of them 
are adjacent, then in the transformed graph 7(P) those points each be- 
long to at least m different uniformly directed cycles that have no com- 
mon sides, where m is the parameter in T; (I, p. 334). The m cycles have 
the 7 points fi, fo, ... , fim common. If some of the points are adjacent, 
then in the transformed graph T(P) each of the points fi, fo, ..., f, again 
belongs to at least m uniformly directed cycles, but those m uniformly 
directed cycles have as many sides in common as there are pairs of adja- 
cent points in the set fi, fo,..., fr 

Proof. Operations T; and 7; result in the fusion of all corresponding spe- 
cialized points of the m primordials (I, pp. 335-36). Hence each point f; 
(¢ = 1,2,..., 7) in T(P) belongs simultaneously to m uniformly directed 
cycles. If a pair of points f;, f, are adjacent, then operations 7; and T3 also 
result in the fusion into one line of all the m lines f; > fx. Hence f; > f;, 
is a common side of m cycles. 

Corollary. If in the primordial graph P a uniformly directed cycle con- 
tains only the specializable points fi, fo, . . . , f,, then those points belong 
in T(P) to only one uniformly directed cycle. 

Proof. The m cycles each consist of r sides. Since there are r adjacent 
pairs of points, therefore the m cycles have r sides in common; in other 
words, they fuse into one cycle. 

The theorem and corollary hold, as is readily seen, for non-uniformly 
directed as well as non-directed cycles also. The reason for our emphasiz- 
ing the uniformly directed cycles is as follows. 

In I we have seen (pp. 326-27) that if in a uniformly directed cycle one 
biological function is destroyed, the remaining biological functions repre- 
sented by the points of the cycle also are destroyed, since each following 
biological function depends on the preceding. If in the primordial organism 
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the r specializable points fi, . . . , f, all belong to only one uniformly di- 
rected cycle, then the destruction of a biological function represented by 
any residual point of the cycle destroys all other biological functions in 
that cycle, in particular, the specializable biological functions fi, fa, . . . , 
f,. In T(P), however, except for the very special case of the corollary, de- 
struction of one point in one of the m cycles breaks only the chain of bio- 
logical functions in this particular cycle. There still remain m — 1 intact 
cycles, to each of which the set of biological functions, represented by the 
points fi, f2,...,/, belongs. Thus in this case those biological functions 
are not affected by the injury. Two adjacent points f; and fi41 are joined 
by only one line. But, in general, two points are joined by m lines. We 
may remove at leastm — 1 residual points between any two specialized 
non-adjacent points f; and fi41, still leaving one way fi, — fr41 intact. 
(We can remove more than m — 1 such points if there is more than one 
residual point between f;, and fi11 in P.) Hence if v denotes the number 
of pairs of adjacent points in fi, fo,..., f,, then we can remove at least 
as many as (r — v)(m — 1) residual points and still leave at least one 
uniformly directed cycle intact. 

Translating this topological result into biological language, we may say 
that the greater m, that is, the more developed and differentiated an or- 
ganism, the greater are its possibilities of withstanding the loss of some 
residual biological functions, the role of which is taken over by other cor- 
responding residual biological functions. 

In the same manner as theorem 1, we can prove generally the following 

Theorem 2. If in the primordial graph P there are p directed ways that 
lead from a specializable point f; to another specializable point f;,, through 
some residual points, then in P(T) there are at least mp directed ways 
leading from f; to fy. 

The biological implication of theorem 2 is the same: a larger number of 
directed ways between f; and f; provides a greater safety factor against 
results of destruction of some of the biological functions which are “‘be- 
tween’ the biological functions represented by /; and f:. 

A similar situation holds for the destruction of a specialized biological 
function represented by one of the points of the set fi, fo, ..., fr. In the 
primordial graph P let the point f;, be preceded by the point fy and fol- 
lowed by f;’. These two points may either be residual points or they may 
be elements of the set fi, fo,..- fr. In this case fy = fia; fr = fra 
According to 7; (I, p. 334) in T(P) the point f, has in general a subsidiary 
point, f:,, which is connected to f; and f;’ in the same way as f; is connect- 
ed to them. Hence the removal of f, from the uniformly directed cycle still 
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leaves the other points fi, fo,..-, fr-1, fist, --» fr belonging to a uni- 
formly directed cycle, the point f, being “by-passed,” and the directed 
way f, >fxr—f, being substituted by fi — fy >f . Thus, as we have 
seen in III, while a higher organism has in general a lesser ability of re- 
generation of a lost biological function, it has a greater ability to adapt 
itself to such a loss, by substitution of other biological functions. This 
again is a well-known fact of biology. 

We may go even further. In I (p. 340) we suggested that when the dif- 
ferentiation is completed, that is, when m = m (where 1 is the total num- 
ber of specializable points), we may apply the transformation T to this 
“completed” graph 7(P) as we applied it to the primordial graph P. In 
this way we describe specializations between various individuals and ob- 
tain a topological representation of social relations. We may, however, 
apply the transformation T to any T(P), not necessarily to the completely 
differentiated one. This corresponds to actual facts, since social differen- 
tiation is found among some relatively lower animals. The specialization 
may be either a social one, a biological one, or both, different individuals 
of the species being biologically different for the performance of different 
biological functions. The most common example is the sexual differentia- 
tion, which in human society is both biological and social. Sociology may 
justly be considered as the ecology of human society. 

Purely biological specializations occur at rather low levels of develop- 
ment. Thus plants specialize in producing organic material from water, 
carbon dioxide, and inorganic salts. Some animals use this organic material 
by feeding on plants. An organism is never completely isolated. A large 
part of its environment is formed by other organisms with which it inter- 
acts biologically. Therefore, the representation of an organism by a graph 
that describes the interactions within the organism is an unrealistic ab- 
straction. We should consider the graph of an organism as a partial graph 
of a more general one, which includes other species with which the or- 
ganism interacts. Since the number of species interacting with a given one 
is finite, though sometimes very large, we still are dealing with finite 
graphs. 

We may now apply the transformation T to this general graph as we 
applied it to the primordial. Since the theorems 1 and 2 of III and the 
theorems above in this paper are independent of the choice of the graph 
to which we apply 7, we therefore arrive at the same results, only with 
a much broader biological scope. The removal of some specialized points 
now may mean the destruction of some of the organisms on which the 
given organism feeds, or of some particular food component. Interpreted 
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biologically, theorems 1 and 2 above again mean that the higher an or- 
ganism, the more adaptable it is to changes in biological environment. 

Thus regardless of the choice of the primordial graph, the transforma- 
tion T leads us to two well-known biological laws: that of decreased re- 
generating ability and that of increased adaptability with increasing de- 
velopment. As we have remarked in I, the transformation T does not de- 
scribe some important biological facts as will be made clear below. There- 
fore, the finding from it of two actual biological laws must be considered 
as coincidental. We may, however, again, as in III, pose the interesting 
question: What conditions must the transformation satisfy in order to 
lead to given known biological laws? The answer to that question will 
narrow the possible choices of a proper transformation. 

Within this given system of abstract biology we may, however, find 
other general laws by further study of the general properties of T. 

Different types of general laws will be obtained by studying such topo- 
logical properties as are invariant with respect to T. (The author owes this 
suggestion to Dr. Ernesto Trucco.) The two laws obtained so far deal with 
differences between organisms of different degrees of development. The 
laws obtained from the study of invariance will emphasize the similarities 
between different organisms. 

The following example illustrates a still different kind of biological rela- 
tion which may be derived in this system of abstract biology. 

It is natural to interpret m as the number of organs of an organism. For 
each of the m component primordial graphs contains some specialized 
points as well as their subsidiary points and forms a natural unit of the 
graph T(P). Each of the x — 1; subsidiary points is attached to a residual 
graph. In the interpretation which we have given in I those residual 
graphs correspond to different tissues of the organ. Therefore the total 
number 7 of different tissues in the organism is equal to the total number 
of residual graphs in 7(P). But that latter is equal to m + n (m — 1) (I, 
p. 335). Hence we have the relation 


tT=m+n(m—1), (1) 


which connects the total number of different tissues with the number m 
of organs and the number » of specializable biological functions. 

So much for the transformation T. If we add to T a definite choice of 
a primordial graph, then, as we have seen, we obtain the full topological 
description of all possible organisms, as well as the number of all possible 
different organisms. 

Thus far our considerations have been of a static nature. A description 
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of the development of the organic world requires, as we pointed out in I 
(pp. 339-40), an ordering of the sequence in which different biological 
functions are specialized. From the topological point of view it would be 
desirable to make this ordering depend on some topological characteris- 
tics of the points that represent the different biological functions. What- 
ever the solution of this problem will utimately be, such an ordering 
will give us the time-sequence of the development of the organic 
world. 

Similar considerations may be applied to individual organisms. Actual- 
ly, the ovum, from which an organism develops, cannot be identified with 
the primordial organism. But in our system of abstract biology we may 
study the situation in which the cell from which a metazoan develops 1s 
the primordial organism. In that case the ordering of the sequence of spe- 
cialization leads to a topological description of embryology and ontogene- 
sis. We may consider a cycle that begins with the primordial graph and 
leads then to a sequence of more and more complicated graphs T(P), 
characterized by larger and larger values of m (<m) and by a definite 
sequence of the other parameters in 7. When the sequence has progressed 
to a point characterized by a given set of parameters, it stops, and the 
process begins again with P, thus symbolically representing the develop- 
ment from ovum to an adult organism, which in its turn produces ova. 
The so-called biogenetic law is then obtained by postulating that the se- 
quences in the philogenetic and ontogenetic developments are determined 
by the same topological criteria. Only an approximate validity of the bio- 
genetic law would follow from environmental differences in the develop- 
ment of the organism and that of the species (see p. 212). 

Criticism of the above system. Suggestion for a different approach. As a 
system of abstract biology, the one discussed above is as good as any 
other. However, inasmuch as our goal is to find a system which is iso- 
morphic with the observed biological world, we must discuss the above 
system from that point of view. Then we find that it does not represent 
the actual situation correctly. 

While in actuality the specialization of cells and tissues frequently hap- 
pens with a loss of some other functions, very much as described by the 
transformation 7, yet there is an all-important type of specialization 
which is of a very different nature. 

The process of ingestion is, for example, specialized already in the 
Coelenterata. Yet in a hydra the ingestion is accomplished both by some 
cells in the gastro-vascular cavity, as well as by the animal as a whole. 
The motions of the tentacles and of the whole body of the hydra produce 
that latter ingestion. In every animal every cell retains the function of 
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elimination of waste materials. Yet, beginning with the Platihelmintes, 
there is a special organ of elimination that does not consist of cells that 
specialize in elimination, since every cell does so. The flame cell, an essen- 
tial part of this organ in Platihelmintes, is a cell which specializes in pro- 
pelling a liquid. Again, every cell of every organism retains the function 
of respiration, yet in higher animals we have special organs like gills or 
lungs which perform the respiration for the organism as a whole. The lungs 
and the respiratory passages are not made of cells which specialize in res- 
piration. The respiration of the organism as a whole is performed by differ- 
ent tissues, some of which specialize in contractility, like the musculature 
which contracts and expands the chest cavity, some in catching particu- 
late matter in the inspired air, like the hairs which line the respiratory pas- 
sages, etc. And so it is with other major special functions, such as diges- 
tion, defecation, and internal transport. In higher animals, digestion and 
defecation functions are lost by all cells and occur only extracellularly. 
Yet these specialized functions of the organism as a whole stand to each 
other in the same relation as do the corresponding non-specialized func- 
tions of each cell. 

If we try to interpret this situation in terms of graphs, we must look 
for such transformations as lead to graphs in which different sets of par- 
tial graphs stand to each other in the same relations as do some points of 
those partial graphs to each other. Such transformations can be found, 
and their study should be our next step, if we hold to the graph as the 
topological representation of an organism. 

Though we were led to the choice of a graph as a topological representa- 
tion of an organism in a rather natural way, we nevertheless must remem- 
ber that within the principle of biotopological mapping other representa- 
tions may well be possible. The construction of graphs with the property 
mentioned above has something artificial and ad hoc about it. Is it not 
possible that by using point set topology rather than combinatorial topolo- 
gy we may find already known topological systems that exhibit the nec- 
essary property? That this may be the case is perhaps suggested by the 
following consideration. In point set topology we may consider sets M; 
having some common properties, and then we may consider a set A of all 
M's. The different M,’s which are the elements of A may stand in the set 
A in the same relation to each other as do some elements of each M;. Such 
situations seem to be more natural in point set topology. Inasmuch as the 
combinatorial and the point set theoretical approaches in topology are 
basically equivalent and represent a difference in methodology, it may 
turn out that the point set topological approach will be heuristically 
more valuable in topological biology. 
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The description of an organism in terms of a graph, its organization 
chart, appears to be the natural thing to do, and the combinatorial ap- 
proach leads us already, as we have seen, to some interesting results. How- 
ever, it is dangerous in the development of a new field of science to follow 
only one possible way, at the exclusion of others. Therefore, if for no other 
reason, we should attempt the use of a point set topological approach and 
describe an organism not in terms of a topological complex but in terms 
of a topological space. Such an approach may turn out to be of particular 
advantage in describing the primordial organism. 

The graph of the primordial organism, as suggested in I, has, apart from 
many other shortcomings, one basic disadvantage: the biological func- 
tions, which are represented by the different points of the graph, are con- 
sidered as given. They are denoted by such biological terms as ingestion, 
digestion, respiration, etc. The use of this terminology implies that the 
organism and its basic relations are considered as given. It would be much 
more satisfactory if we could describe the basic biological functions and 
their relations in terms of some well-known topological concepts, just as 
some basic concepts of physics are described in terms of the metric of the 
four-dimensional universe. In other words, we should look for known rela- 
tions in topology which are isomorphic with known relations in biology. If 
we can do that, then by a proper renaming, that is, by a proper substitu- 
tion of biological terms for topological, we shall obtain a correct descrip- 
tion of different biological relations. 

What is the logical character of the basic biological relations? And 
which of those relations are basic? A careful consideration shows that the 
essential property of an organism, which distinguishes it from the inor- 
ganic world, is the following: an organism selects from its environment the 
proper parts of this environment which are needed for the building-up of 
the organism and eventual reproduction. It does select the proper parts 
from very different environments, where those proper parts are in differ- 
ent combinations with other parts that are not needed by the organism. 
After the selection, the organism assimilates those selected parts, incor- 
porating them into itself. But the organism plus those incorporated parts of 
the environment remains basically in the same relation toward the envi- 
ronment as the organism was originally. The organism plus those parts 
continues the selection and assimilation. 

In the simplest organism, like bacteria, this selection occurs on the 
molecular, biochemical level. In higher organisms it occurs also on differ- 
ent macroscopical levels. Basically the stimulus-response systems of any 
organism are selector mechanisms. Directly or indirectly we use our 
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senses to select our food or to avoid dangerous situations; the latteris also 
basically a selection act. 

In addition to the selection and assimilation of parts of the environment, 
processes of breakdown occur in every organism, providing, directly or 
indirectly, the energy needed for selection and assimilation. 

In the totality of different elements or units which constitute our uni- 
verse, the organism is a subset of this totality of elements. (Element is not 
used here in a chemical sense, but rather in a point set theoretical one.) 
What does the act of selection logically imply? It implies a division of the 
elements which constitute the environment of the organism in two classes: 
the needed ones and those to be rejected. Selection logically implies the 
existence of rejection. We may thus say that the presence of an organism 
induces in the environment a division of all the elements of this environ- 
ment in two classes: class A of elements that are assimilable, and class B 
of those rejected. The addition of a subset A’ of some element of class A 
to those of the organism M leaves M + A’ in the same relation to the en- 
vironment as is M alone. We shall denote this property as “Property I.” 

Cases where subsets of topological spaces induce a division of the 
points of the space in two classes are well known. The simplest example 
is that of a subset of the one-dimensional space of real numbers, when that 
subset consists of only one point. It induces the Dedekind cut. A Jordan 
curve considered as a subset of the Euclidean plane E? induces the divi- 
sion of the points of that plane in two classes: the inner and outer points. 
Those two examples do not possess, however, the above-mentioned Prop- 
erty I. The following third example is somewhat nearer to what we need: 

In a space S consider a subset M which is not closed relatively to S. It 
induces the division of all points S — M in two classes: those that are 
limit points of M (class A), and those that are not (class B). Since a sub- 
set remains not closed as long as not all of its limit points belong to it, 
therefore, adding any proper subset A’ of A to M still leaves M+ A’.a 
not closed subset, which induces a division of S — (M + A’) in two 
classes: those that are limit points of M + A’ and those that are not. We 
have here an isomorphism between a biological set of relations and a 
purely topological one. 

Under some circumstances M@ may be not connected. If, for example, S 
is the one-dimensional space of real number and / is the subset of all ir- 
rational numbers in the interval (0, 1), then M is both not closed and not 
connected. The addition to M of any finite subset A’ of the set of rational 
numbers in the interval (0, 1) still leaves M both not closed and not con- 
nected. Being not closed and not connected are two closely related prop- 
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erties of M in this case. We have here an analogy with the above-men- 
tioned biological fact that a break up of the organism is biologically an 
essential part of the processes of selection and assimilation. 

This example suggests a possible way for an application of point set 
topology to biology. This way may well lead to a geometrization of biolo- 
gy, similar to the geometrization of physics. Different modalities (the 
word modality is used in the psychobiological sense) of stimulus-response 
relations may perhaps be represented in multi-dimensional spaces, each 
modality corresponding to a dimension. Once a description of a primordial 
organism in point set topological terms is achieved, a proper transforma- 
tion, which produces more complicated topological spaces that can be 
continuously mapped onto the topological space of the primordial, will 
give us a description of other organisms. 

When it was realized that phenomena of the physical universe are iso- 
morphic to some metric properties of a particular abstract four-dimen- 
sional geometry, the theory of the abstract geometries was well developed, 
and the results of that theory could be readily taken over by the physicist. 
From a mathematical point of view it would perhaps seem more aesthetic 
to find a proper representation of life by topological systems which have 
already been studied by mathematicians, due to the purely mathematical 
beauty of such systems. On the other hand, Gauss is supposed to have 
said that he received many an inspiration for his purely mathematical 
discoveries from the study of problems of physics. Also in the case of rela- 
tivity theory, though, as remarked above, the basic geometric theory was 
available to the physicist, the further requirements of the theory led to 
purely mathematical generalizations, such as the introduction by Einstein 
of an antisymmetric fundamental tensor for the possible interpretation 
of the electromagnetic field. 

Whether topologists will provide the mathematical biologists with a 
proper abstract system, ready to be applied, or whether the study of life 


will inspire topologists to further discoveries in their science, only the 
future will tell. 
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A refinement of a model of the vascular system (Cohn, 1954) is discussed. We now consider 
that the major vessels arise from the aorta as secondary branching vessels rather than by 
successive dichotomies of the aorta as previously considered. It is shown that this change does 
not affect calculated values of flow, the number of branchings, or the relative radii of vessels 
at a branching. A comparison is made of theoretical predictions with experimental data and 
good agreement is found. 


The method of analysis of this paper is used to extend the treatments of the previous paper 
to a more general case. It is shown that the results previously obtained are valid. 


In a previous paper (Cohn, Joc. cit.) a model of the mammalian vascular 
system based on engineering principles was proposed. The model was 
based on considerations of a good way to supply a volume of space with 
a uniform vascular system. Basically it was desired to start with a single 
vessel of supply (the aorta) and by successive branchings ultimately sup- 
ply a fine network of vessels (the capillaries) uniformly distributed 
throughout the space. A limit was placed on aortic radius based on tur- 
bulence conditions* and the density of capillaries was limited by oxygen 
diffusion requirements. The intermediate branching system was consid- 

* We may note here a point not made clear in the first paper (Cohn, Joc. cit.). When the 
aortic radius has been limited by turbulence considerations, turbulence will not arise in any 


part of the branched system. This will be made clear in the following. 
The Reynolds number is defined [cf. Cohn, Joc. cit., eqs. (1) and (2)] as 


in which 7 is the fluid viscosity, o is its density, C is the flow through the vessel, and r is the 
radius of the vessel. Then at a point of equal branching at which the radii of the branch ves- 
sels are .794 times the radius of the main vessel (the branched system we have considered) we 
find that the Reynolds number of the branched vessels is 
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It may be seen that the Reynolds number after the branching is always smaller than before 
the branching. Thus, if the largest vessel of the system (the aorta) does not have turbulent 
flow, none of the smaller vessels will. 
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ered as generated by successive equal branchings of the main vessel 
(aorta) and the vessels generated by it. The resistance of the model system 
was then calculated. Good agreement between the model system parame- 
ters and anatomic parameters tabulated by H. D. Green (1950) was 
found. 

One difference between the proposed model and the anatomical situa- 
tion is that in the animal the major vessels arise from secondary branch- 
ing rather than from successive equal branchings of the aorta. Similarly 
the vessels these generate also arise as secondaries from the main vessel. 
This type of vessel branching would originate from an elongated model 
rather than the compact cubical form we originally considered as an ab- 
stract model. The mammalian form may be more accurately represented 


Ficure 1 


by a cylindroid trunk with cylindroid appendages originating from the 
trunk. To supply blood to such an extended form it is reasonable to start 
with a main vessel running the length of the organism with secondaries 
arising from it to supply it various parts. 

With these considerations our first crude model of the vascular tree 
would be replaced by the following. The larger vasculature arises as sec- 
ondary vessels from the aorta. This type of branching eventually is re- 
placed in the smaller vasculature by successive equal splitting of the ves- 
sels. Thus the refinement of the model considered in this paper is only ef- 
fective from the aorta to the medium-sized vessels. From the medium 
sized vessels to the smallest vessels, branching by equal splitting is con- 
sidered. 

As a model of the branched system for the larger vasculature we will 
consider the vessel system of Figure 1. We let AB be the main vessel with 
m secondary vessels of radii 7; and lengths 7; branching from it. If AB is 
the aorta, these secondaries would be the vessels supplying the arms, legs, 
etc. Assume the secondaries give rise to vessel systems of equal resistance 
to flow. This may be justified in the model under consideration by the 
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fact that the aortic resistance is very small compared with the resistance 
of the smaller vasculature. Thus the flow through the secondaries is deter- 
mined by the resistances of the smaller vessels. Since the flow is not deter- 
mined by the aortic resistance, and does not vary as we change the aortic 
parameters, it is justifiable to assume equal flow through all secondaries. 
The anatomical situation justifies this simplication, for the main vessels 
leading from the aorta carry approximately equal quantities of blood. 
We may now evaluate the radii of the system by a method similar to but 
more rigorous than that used previously (Cohn, 7bid., p. 69). 

It is desired that the resistance of the system to flow be low so that the 
heart will do a minimum of work in order to circulate the blood. This low 
resistance could be accomplished by indefinitely enlarging the radii of the 
vessels, but this would be very uneconomical, for the body is limited in size. 
Thus we place a limit on the amount of material (tissue) that may be used 
in construction of the vessel system and work with that restriction. To 
determine the amount of material used we must know the vessel wall 
thickness. As done previously (Cohn, zbid., p. 69) we assume the vessel 
wall thickness to be proportional to the radius. Then 


6;= a7; and A;= ak; , (1) 
in which 6,, A; are the wall thicknesses of the 7th secondary and main 
vessels respectively, 7;, R; are their radii, and a is a constant of propor- 
tionality. 

We assume a constant amount of tissue available for the vessel walls 
of the system. Thus we obtain as a limiting condition on the system that 


M= 27RyLiA1 + 2maRoloAot+. 5 t+ 2arihoit 20 Volo62+. dle 95 


in which M is the total mass of the system. 
Substituting (1) in the above we obtain 


M= IreR?L, + 2roRiL, +. boxer + lrarht... . 


The last equation may be written: 
MO ERAT 2 
an 2 (RLi+ rils) , (2) 


in which the summation is over the 2 secondaries arising from the main 
vessel. 

The total energy loss AE due to fluid flowing through the system may 
be expressed by means of the following equation 


AE => (FiQ;4+ fiws) 5 (3) 
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in which F; and f; are the flows through the ith main vessel segment and 
secondary vessel respectively and Q; and w; are their resistances to flow. 

We must next determine the flows through the various vessels and their 
resistances to flow. We have already assumed that the main vessel gives 
rise to-# secondaries. Thus if /; is the total flow entering the main vessel 
at point A, each secondary carries an outflow of (1/)F,. Then 


fe=—Fi. (4) 


From this it follows that the flow carried by the ith aortic section is the 
initial inflow F, decreased by all previous outflows. This is represented by 
the following equation 


; 1 
Ei ie G-1) 7h 


or 


F; ae “ Fy ° (5) 
To obtain the resistances to flow of the various vessels we use the 
Poiseuille law and find 


Oy and. ee (6) 


in which 7 is the fluid viscosity and o its density. 
Now we may substitute flow and resistance values in our energy loss 
equation. Introducing (4), (5), and (6) into (3) we obtain 


iin he > (2 Ga pce ete 


R‘ Tron? rs 


(7) 
= ie ah 


The constant mass of the system may be introduced by solving (2) for 
R, from which we obtain 


MSN Pt Sn 
ae DRL: ay (8) 


tak 
| 


L; 
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Introducing (8) into (7) we obtain 


2 273 
TON = n Ly 
SaF | M 


27a 


ta a >> [ (1a it] - (9) 


n 
3 p2 2 =e ; 
—2 L:—274,| a w=2 u 0 


In order to obtain the relative radii of the system we take partial deriva- 
tives of AE with respect to R; and r;. This differentiation leads to 


DAE | oe occas em (10) 
OR, : x To Nn? 
and 

dAE [nn 178nFi 

alee a 


We may obtain the radii for minimum energy dissipation by setting 
the above equations equal to zero. This gives 


From (10), (11), and (12) we then obtain 
i= [1-4], (13) 
rm. (14) 


We shall now interpret these results and relate them to our previous 
model of the vascular tree. In the mammal about eight major vessels origi- 
nate from the aorta, which courses along the length of the trunk. In gen- 
eral, each of the secondaries also courses along a cylindrical volume, limb, 
or gut, etc. Thus the secondaries branch by giving rise to their own sec- 
ondary vessels. Finally, however, a stage is reached where the vessel must 
supply a compact volume of tissue, and then we may use our former 
method of analysis involving equal vessel branchings. 

We have previously shown (Cohn, ibid., p. 65) that a single vessel gives 
rise to eight vessels after three equal branchings and that at each split 
the radii of the daughter vessels are .794 times the radius of the previous 
vessel (Cohn, ibid., p. 70). Thus after three splits the radii of the vessels 
are (.794)3 = .500 times the radius of the original vessel. Substituting 
m = 8 in (14) we find that our present calculation gives the same radius 
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for each of the eight secondary vessels arising from one main vessel as 
would be obtained with equal successive branchings after three branch- 
ings. Thus the main factor affecting flow in the system (the vessel radius) 
is the same whether we consider branching by equal splitting or by sec- 
ondaries arising from one main vessel. 

We emphasize again that the manner of branching considered in this 
paper takes place for the first two, and at most three, generations of ves- 
sels. Since the radii are the same in either method after this number of 
branchings, and since the main component of resistance to flow arises from 
the last generations of branching (the arteriolar and capillary level) we 
find that the calculations of flow values and number of vessels made in 
the previous paper (Cohn, ibid.) are not significantly affected by the new 
type of branching here considered for the first two or three generations of 
vessels. 

Before summarizing the results obtained so far we may note that the 
conditions by which the relative radii of the system have been obtained 
are of a more general nature than those of the preceding paper (Cohn, 
ibid.). The present treatment also corrects an error of the previous paper. 
There we considered the mass as being constant for a single branching re- 
gion, a situation which is physically meaningless. In the present paper the 
mass of the entire system has been limited and the energy dissipation of 
the entire system considered. The radii of the entire system have then 
been determined by minimizing the energy dissipation of the system. This 
method of analysis may be applied to the problem treated in the first 
paper without altering the results obtained in that paper. In the following 
we shall apply the present corrected treatment to the case of successive 
equal branchings which was studied in the previous paper. 

Consider the system illustrated in Figure 2. The resistance © to fluid 
flow of a single vessel is given by the Poiseuille relation [cf. eq. (6)] 

8n 1 


To r*’ 


Q (15) 
in which 7 is the viscosity of the fluid, its density, J the length of vessel, 
and r the radius of the vessel. 

The total resistance of the system Qr is 


Qpr =A + FM +404... 0, (16) 


where Q; is the resistance to flow of a single vessel of the ith generation. 
The factors before the Q; values are the inverses of the number of vessels 
of that generation. 
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Substituting from (15) in (16) we obtain 


39, 18a k iL, Wap ths 
ae Begg Tica Ma se pit = D8 ee 7a! (17) 
which may be written 
87 Lo 8 . 1 l; 
—— 
Ug To ries anes (18) 


Before minimizing the resistance to flow we must consider the total 


FIGURE 2 


mass, M, of the system which may be obtained as in the previous paper. 
The expression for total mass is 


M=2raro+2-Irarih+...2"+2raraly . (19) 
From this we obtain the following expression for 75: 
ro 2*r; 20 
ie ib ee ss rid J. EN) 
Substituting from (20) into (18) we obtain 


87 i 8n~x 1h 
(21) 
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We may now differentiate this expression and obtain 


OQ7 _ 8 I geass 8 l; 
ir | gle De wea (22) 


Setting this partial derivative equal to zero (the condition of minimum 
resistance) and simplifying we obtain 


re= 2G: (23) 
This result is a more general derivation of the result of the first paper in 


which it was found that the daughter vessels at a single branching have 
.794 times the radius of the vessel from which they originate. The simi- 


TABLE I 
PARAMETER THEORY EXPERIMENT 
Aortic radius .43 cm. 38) O00 
Diameter of vessel after split rela-| .794 about .8 
tive to original diameter 
Total blood flow! 57.5 gm..sec.| 40 gm. sec. — 
Number of classifications of arteries? | 11 9 


1. We considered the total flow as given in determining the aortic radius by turbu- 
lence considerations (Cohn, ibid., p. 71). However, with the aortic radius so obtained we 
constructed the remainder of the vascular tree. Then, by imposing the known arterial 
pressure of 100 mm. Hg. we obtained the flow figure listed above. The fact that this figure 
is close to that used to obtain the aortic radius shows that we have made a reasonable as- 
sumption for radius size, but the flow figure given above is independent of the original 
value used. In a similar manner we might have assumed the flow value experimentally 
ee (40 gm. sec.~1) and from this derived the blood pressure necessary to achieve this 

ow. 


2. The discrepancy in listing of major types of arteries is mainly a matter of classifica- 
tion. Green lists 40 major arteries, as arising from the aorta. From this it is obvious that 
he is considering not only those directly connecting with the aorta, but also some arising 
as branchings of these. Thus his classification is somewhat broader than ours, and so the 
number of classifications he considers is smaller. 


larity of the two results may be noted if we call 7;,, the daughter vessel 
radii and 7; the radius of the vessel from which they originate. Then 


Me eC aa 


tS A238 Tijats 
which reduces to 


“t= 794, (24) 
which is the result previously obtained. 

Since several questions concerning our original model have been re- 
solved, this is an appropriate place to summarize the results obtained so 
far. As done before (Cohn, ibid.), we will compare our results with anatomi- 


cal data for a 13 kg. dog as assembled by Green (ibid.). This comparison 
is made in Table I. 
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We may note the good agreement between theoretical and observed 
parameters of the arterial system. On minimum energy dissipation con- 
siderations we have constructed a vascular system very similar to the 
anatomical system. This is interesting in that it illustrates one principle of 
development of the vascular system that probably could be applied to the 
analysis of other biological systems. We have also obtained analytical ex- 
pressions for the dimensions of all component parts of the vascular system 
which might be of interest in the solution of flow problems, reflections of 
pressure waves, etc. Finally, if the present results warrant further refine- 
ment of this model it might be of interest to consider secondaries arising 
from the aorta carrying unequal flows, or the effect of a small vessel aris- 
ing as a secondary among larger secondaries. 


The author would like to thank Professors Rashevsky and Landahl for 
their suggested extensions of the previous paper which led to the develop- 
ment of this paper. 

This work was aided in part by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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The information content of an organism determines to a large extent its ability to perform 
the basic vital functions: selection of food, breaking up of the food molecules into appropriate 
parts, selection of those parts, and their assimilation. The information content needed is very 
large and requires a sufficiently large complexity of the organism. The information content of 
an organism is largely determined by the information content of the constituent organic mole- 
cules. The information content of the latter is in its turn determined by the number of physical- 
ly distinguishable atoms or radicals of which the molecule is composed. The different arrange- 
ments of atoms in a molecule are represented by the structural formula, which is basically a 
graph. It is shown that the topology of this graph also determines to a large extent the infor- 
mation content. Different points of a graph may be physically indistinguishable; in general, 
however, they are different in regard to their topological properties. A study of the relations 
between the topological properties of graphs and their information content is suggested, and 
several theorems are demonstrated. A relation between topology and living processes is thus 
found also on the molecular level. 


Information-theoretical problems in biology have been recently raised 
by a number of authors (Quastler, 1953; Morowitz, 1955). The following 
information-theoretical aspect of very basic biological phenomena, though 
not explicitly stated, seems to be implied in the interesting considerations 
of J. von Neumann (1951) on self-reproducing automata. 

Perhaps the most basic characteristic of living organisms is their ability 
to construct other similar or even more complex organisms through as- 
similation of certain component parts of their non-living environment. 
The organism, whether multicellular or unicellular, or even a hypothetical 
“living molecule,”’ does not have available in the environment its com- 
ponent parts in a “pure form.” Those component parts are contained in 
the environment as parts of other non-living entities. The organisms must 
break up those non-living entities into smaller units in a proper way, and 
then again select from this debris the component units needed for building 
other similar living organisms. After such a selection the organism ar- 
ranges the component units in a proper manner. This is the essence of 
digestion, ingestion, and assimilation—the basic properties of the living. 

Viewed as an automaton, the organism must have, in von Neumann’s 
terminology (Joc. cit., pp. 29-31), instruction for the proper break-up of 
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environmental units, the proper selection of the necessary components, 
and their proper arrangement. From the point of view of information 
theory, the organism must have a sufficiently large information content. 
It must contain all the information necessary for the proper break-up of 
the environmental unit (digestion), the selection of proper units resulting 
from the break-up, and the building of a duplicate organism from those 
units (absorption and assimilation). 

Leaving aside the different, sometimes divergent, estimates of the infor- 
mation contents of cells and organisms (Quastler, loc. cit., pp. 251, 263), 
we can see offhand that the information content must be very large. In 
order to have a sufficiently large information content, the organism must 
be sufficiently complex. Therefore the above considerations put a lower 
limit on the complexity of even the simplest organisms. 

This complexity may be to some extent measured by the number of 
different, distinguishable units of which the organism is composed. If we 
have any large number of indistinguishable units, the information content 
of each unit and of the whole aggregate is zero. If we have m distinguish- 
able units with equal probability of choice, then the information content 
is equal to loge 1. 

If we consider as elementary units the different kinds of atoms of which 
living organisms are made, we find that this number is very limited—of 
the order of 10. Therefore second-order units must be made from combina- 
tions of those elementary units. Such second-order units are the diverse 
kinds of molecules of which an organism is composed. The situation is 
quite analogous to the formation of a very large number of distinct words 
from a limited number of letters in the alphabet of a written language. 

From the point of view of this analogy the above problem of minimum 
complexity becomes formally analogous to the following one: 

Given a written language, how large, in terms of the total number of 
words, must a book printed in that language be in order to contain com- 
plete information necessary to manufacture this book? 

The answer depends on the nature of the environment, on what is given 
for the manufacturing. If a complete linotype machine, bookbinding 
machine, paper, ink, etc., are available, the book need not be too large. If, 
however, only rough materials are available out of which one must make 
the paper ink, machines, ete,, then the size of the book will be beyond 
the practical limits of a single volume, and we shall have the problem of 
reproducing not a book but a whole library. In the case of the living the 


situation is nearer to that of the second alternative; hence the tremendous 
complexity of organisms. 
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The analogy between the information content of a written text and 
that of some organic molecules has been pointed out by H. Branson 
(Quastler, Joc. cit., pp. 84-118) and more recently by G. Gamow (1954a, 
b). The physically different atoms, or their physically different aggre- 
gates, play the role of letters and words. 

It is, however, possible for an aggregate of identical, physically indis- 
tinguishable units to have a large information content. The units in this 
case, though indistinguishable physically, are different through the differ- 
ence of their relations to each other. 

A structural formula of organic chemistry is basically a one-dimensional 
topological complex or graph, whose points are occupied, in general, by 
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physically different and distinguishable atoms. The theory of graphs has 
been successfully applied to some problems of organic chemistry (Kénig, 
1936). 

However, even if the points of a graph are physically identical and in- 
distinguishable, the graph may possess a finite, possibly large, information 
content. The following examples illustrate this. 

In the non-directed graph of Figure 1a, there is no way of distinguishing 
between the points 1 and 2. The information content of that graph 1s zero. 
In the directed graph, 1b, the two points 1 and 2, though physically in- 
distinguishable, are topologically different: point 1 is an end point of a 
directed line, point 2 is a beginning point. The information content of this 


graph is, therefore, 1 bit. 
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In the undirected graph, 1c, points 1 and 3 are indistinguishable, both 
being of degree one. Point 2, however, is different, being of degree two. 
The probability of a point in 1c chosen at random being of the first type 
is 2, the probability of it being of the second type is 3. Hence the informa- 
tion content of graph Ic is 


I,= — (2 loge +3 log, 3) = 0.92. (1) 


If, however, from the graph 1c we make the directed graph 1d, then all 
three points become topologically distinguishable: point 1 is of degree 
one and is an end point of a directed line. Point 2 is of degree two and is 
an initial point of one directed line and an end point of another directed 
line. Point 3 is of degree one and is an initial point. The information con- 
tent of 1d is log, 3 = 1.59. 

In both the undirected and the directed regular graphs 1e and 1f all the 
points are indistinguishable, and their information content is zero. In the 
graph lg all three points are topologically different, and its information 
content is logs 3. 

Now consider graph th. It has 3 points of degree two (1, 3, and 4), and 
two points (2 and 5) of degree three. However, points 1 can be distin- 
guished topologically from either point 3 or 4, while the latter two are in- 
distinguishable. Point 1 is adjacent to two points of degree three, while 
each of the points 3 and 4 is adjacent to one point of degree three and one 
of degree two. Hence we have 3 distinct classes of points, 1; 2 and 5; and 
3 and 4. The corresponding probabilities are 4, 2, and 2. Hence the in- 
formation content of ih is 


Let us, however, add an extra point to graph 1h so as to obtain graph 1i. 
Now not only points 2 and 5, but also the two points 3 and 4, which are 
indistinguishable in 1h, become different. For although both points 3 
and 4 are of degree two, yet one of them (point 3) is adjacent to one point 
of degree four and to one of degree two, while the other (point 4) is 
adjacent to one point of degree three and one of degree two. Moreover, 
point 3 is two units distant from point 6 of degree one, while point 4 is 
three units distant from 6. The information content of 1i is 


I;=log, 6=2.6. (3) 


Adding to 1i another point, so as to obtain graph 1k, we find seven - topo- 
logically distinct points, and an information content 


I, = log, 7 = 2.82. (4) 
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On the other hand, adding to 1i an extra point so as to obtain graph 11, 
we reduce the information content, because now the points 2 and 5 are 
again indistinguishable. So are points 6 and 7, as well as 3 and 4. The in- 
formation content is now 


I, =+4 log, 7+$ log = 1.96. (5) 


The information content of all the above graphs is determined by purely 
topological properties. We may call it the topological information content. 
If, in addition, some of the points are physically distinguishable, then the 
total information content is composed of the topological information con- 
tent and of the information content due to the physical differences. Thus 
if in graph 1b we make the two points 1 and 2 physically distinct, as in 
Figure 2, the information content becomes log, 4 = 2 bits. For we have 


e——<—9 
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four equally probable possibilities: white initial point, white end point, 
black initial point, black end point. 

The above considerations raise rather interesting purely mathematical 
problems. What general relations exist between the topological properties 
of a graph and its topological information content? How many different 
graphs of points each are there with the same topological information 
content? The study of such problems is likely to lead to a number of 
theorems which may be both of purely mathematical interest and of value 
in organic chemistry and biology. 

To illustrate the kind of mathematical problems that such a study will 
involve, we prove here a few simple theorems. 

Theorem 1. The topological information content of a finite connected 
regular non-directed graph is zero. 

Proof. A finite regular connected non-directed graph is either a polygon 
or is homeomorph to a graph formed by the edges of one of the five regular 
polyhedra. The identity of all points is readily seen. 

Theorem 2. If to a non-directed graph which consists of a cycle with 
2n + 1 points and has, according to theorem 1, a zero topological informa- 
tion content, we add a point, joining it to one of the 2m + 1 points, the 
new graph thus obtained has a topological information content equal to 

1 


n 


net 108 + 1) = (6) 
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Proof. The new graph has 2n + 2 points. One of them, the added one, 
is of degree one; the point to which it is joined is of degree three. The re- 
maining 2n points are divided into x pairs of identical points, the points of 
each pair differing from those of other pairs by their topological distances 
fronr the added point. 

The probabilities of the points of degrees one and three are each 
1/(2n + 2); the probabilities of the other m are each 


Te en 
Qn+2 nt+1° 


From this (6) follows. 

Theorem 3. If to a non-directed graph which consists of a cycle with 2” 
points and has, according to theorem 1, a zero topological information con- 
tent, we add a point, joining it to one of the z points, the topological con- 
tent of the new graph thus obtained is equal to 


2 (m—1) ee 2n+1 


In+1 ay) 


logs (2n+1) + (7) 


2n “ Qn+1 

Proof. The new graph has 2% + 1 points. The added point is of degree 
one; the one to which it is connected is of degree three. The remaining 
2n — 1 points contain 2 — 1 pairs of identical points, the points of 
each pair differing from those of others by their topological distance from 
the added point. Finally there remains one point which differs from all 
others by its topological distance from the added point. The correspond- 
ing probabilities of each of the distinct classes are 1/(2m + 1) for each of 
the three single points and 2/(2m + 1) for each of the m — 1 pairs. From 
this (7) follows. 

Theorem 4. A graph which is a uniformly directed cycle (Rashevsky, 
1955c) has a zero topological information content. 

Proof. All points are indistinguishable, each being an initial point of a 
directed line and an end point of another directed line. 

Theorem 5. If to a directed graph which consists of a uniformly directed 
cycle with points and has therefore, according to theorem 4, a zero topo- 
logical information content, we add a point, joining it by a directed line to 
one of the ” points, the topological information content of the new graph 
is logs (7 + 1). 

Proof. The new graph has 2 + 1 points. The added point is of degree 
one; the point to which it is joined is of degree three; the remaining 
m — 1 points, though all of degree two, are nevertheless all distinct, be- 
cause they all differ in their distances from the point of degree three, 
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counted in the direction of the arrows. Hence there are » + 1 distinct 
points, all of equal probability. 

Connections of a different kind between information theory and 
topology already appear in the work of A. Bavelas (1950) as well as in the 
classical book of C. Shannon and W. Weaver (1949). The importance of 


topology in biology has been shown in several recent papers by the au- 
thor (Rashevsky, 1954; 1955a, b, c). 

The information content of an organism is largely responsible for the 
performance of the basic vital functions: selection and digestion of food, 
absorption, assimilation, and reproduction. The information content of an 
organism in its turn largely depends on the information content of the 
constituent organic molecules. And that information content depends, as 
we have seen, on the topology of the molecule. Thus we again find, this 
time on a molecular level, a close relation between topology and life. 


This work was aided by a grant from the Dr. Wallace C. and Clara A. 
Abbott Memorial Fund of the University of Chicago. 
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BOOK REVIEW 


Pavut F. Lazarsretp, Epitor. Mathematical Thinking in the Social Sciences. 1954. 444 pp- 
Glencoe, Ill.: The Free Press. 


This book is composed of eight articles. All of them deal with interrelations between mathe- 
matics and the social sciences, but otherwise they differ in content and scope. The “building 
of models”’* and the importance of probabilistic concepts for analyzing the field of behavioral 
science can nevertheless be singled out as the two main topics discussed in the papers, some of 
which were originally presented as lectures at a meeting sponsored by the Department of 
Sociology at Columbia University. 

The book under review is not, therefore, a textbook on “Mathematics Applied to the Social 
Sciences” nor even a general survey of what has been done in this field. To quote from Lazars- 
feld’s introduction: ‘The general idea was to represent a number of ‘cases’ as examples for 
situations as they might typically develop between the social sciences and mathematics,’’ and 
again: “.... this symposium makes no claim to having covered in any exhaustive way the 
role of mathematical thinking in the social sciences. It will have served its purpose if its read- 
ers get a feeling for the present state of affairs, and if some of them become stimulated to con- 
tribute their own efforts toward a further development.” y 

The applied mathematician will be interested in this book which points to new possibilities 
for using mathematical tools and ideas. However, the group for which the volume was pri- 
marily intended consists of ‘‘social scientists with only a modicum of training in mathematics.”’ 
Consequently, mathematical complications have been avoided as far as possible, always striv- 
ing to formulate the underlying ideas with the simplest equations. Even so, a mathematically 
untrained reader will undoubtedly have some difficulty in understanding all the details or in 
evaluating the importance and meaning of the indicated calculations, but he should certainly 
be able to follow the basic ideas. Thus, for instance, the reader is helped in comprehending the 
rather formidable looking formulae on pp. 38-39 (in Anderson’s paper) by having the main 
steps constantly verified on a simple (“2 opinions only’’) example. Again, in the fifth paper, the 
detailed mathematical theory of Guttman’s principal components is not given but numerous 
graphs and explanatory remarks are inserted to help visualize the important qualitative fea- 
tures of principal components. Finally, in the last article (by Simon) some results are stated 
without proof; they are based on generally well-known properties of linear differential equa- 
tions with constant coefficients. 

We are thus confronted with eight individual attempts either to handle mathematically a 
problem taken from the general field of social science or to clarify some basic questions con- 
nected with the application of mathematical methods to that field. In what follows we list the 
titles of the eight contributions with a brief summary (if possible) and a few remarks on each. 
The latter are by no means exhaustive, nor should the length of our comments ona given article 
be taken as an indication of that article’s value. 

I. “Probability Models for Analyzing Time Changes in Attitudes,” T. W. Anderson. 

In this paper an “attitude” on a given subject is considered to be a choice from a specified 
set of responses to a certain question, each response being denoted as an “opinion.”’ An indi- 
vidual’s attitude is ascertained at regular time intervals 4, f2, fs, ... . It is assumed that the 


* Cf. also the interesting remarks by Murray G. Murphey in The Saturday Review, Febru- 
ary 5, 1955, page 21. 
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probability of changing opinion at time f, should depend only on the opinion held at ¢, 1. We 
are thus naturally led to the formalism of Markov chains to describe mathematically the 
change of attitudes with time. To extend the model for the case of more than one individual, 
the following two (not very realistic) assumptions are introduced: 

a) all the individuals have the same matrix P of transition probabilities, which, moreover, 
is assumed not to vary with time, and 

b) the individuals change their opinions independently of one another. 

The conclusions reached are interesting from the purely mathematical point of view (being an 
extension of the usual Markov chain theory). As Anderson himself points out, this is a typical 
example of a model constructed “ad hoc.’’ Certain assumptions are translated into mathe- 
matical language and conclusions reached from them by mathematical manipulations. In 
general the conclusions thus obtained will be found to hold true only if the initial assumptions 
were fulfilled. The example given by Anderson in his third section shows this very clearly: for 
the particular case studied it is found that one constant matrix P is inadequate to describe the 
situation. However—and this is a general feature of such models—even discrepancies between 
predictions and empirical results may be helpful in pointing the way toward better theories. 
In Rashevsky’s words (page 68, following article): “. . . the mathematical analysis does in 
many cases enable us to determine how our ‘if? should be modified in order to give the quantita- 
tively correct ‘then.’ ’’? Anderson also indicates possible and plausible extensions of his model. 

II. “Two Models: Imitative Behavior and Distribution of Status,’ Nicolas Rashevsky. 

a) Being confronted with the choice between two mutually exclusive reactions, an indi- 
vidual will be stimulated to perform one of them predominantly if he sees other individuals of 
his group behave in the same manner. It follows from the mathematical treatment that under 
certain conditions (i.e., if the imitative interaction is sufficiently strong) practically the whole 
group will perform one of the two reactions, whereas in the absence of imitation the population 
would have been equally divided in its tendencies toward the two modes of behavior. Some 
interesting results are obtained (e.g., the minimum size of a group necessary for “crowd be- 
havior’’). In contrast to Anderson’s, this model is based on Rashevsky’s own theory of neural 
excitation, but this point is not dwelt upon very extensively. 

b) The second of these lectures deals for the most part with a very general statistical theory 
describing the distribution of social status among individuals of a population and the formation 
of social classes. 

Rashevsky’s considerations about the actual measurement of social status (at the beginning 
of the second lecture), and, more generally, about his studies on neural excitation as a theo- 
retical foundation for psychometric measurements (first lecture, pp. 77-78) should be of par- 
ticular interest to social scientists. 

III. “An Expository Analysis of Some of Rashevsky’s Social Behavior Models,” James S. 
Coleman, 

This paper by Coleman is far more than a simple exposition of someone else’s work. The 
structure and content of the models, the implications, limitations, and meanings of the param- 
eters, the time dependence of certain quantities, the difference between exogenous and endog- 
enous variables, and, finally, the mutual interrelations between the various quantities being 
put together to form the complete model are carefully analyzed according to a well-laid-out 
plan. At the end Coleman shows by means of two examples how a different inteprertation of the 
same equations can be used to account for phenomena totally different from those for which 
the equations were originally intended. 

Coleman’s paper constitutes a real contribution to theoretical social science (at least to its 
“model building” aspects). It will also greatly assist in understanding the immediately pre- 
ceeding article by Rashevsky, which some readers may find hard to follow. 

IV. “Probability in the Social Sciences,’ Jacob Marschak, 3 lectures. 

In his first lecture (entitled “Probabilities and Norms of Behavior’’) the author discusses in 
particular the notion of subjective probabilities which is of special importance to the student 
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of behavioral OTe oe has been the subject of much speculation by mathematicians and 
aad People eee ates and behavior depend on their subjective probabilities (1.e., their 

egrees of belief’’ that some events may happen) and, conversely, the latter can be deduced 
from observed actual behavior, provided an individual acts rationally or “‘as he should.’”’ The 
lecture ends with a simplified presentation of Bayes’ Theorem. 

The second lecture (“Probabilities and Descriptive Social Science’’) emphasizes the neces- 
sity of using statistical methods in the social sciences after having set up some theoretical hy- 
pothesis which we want to test. This is because mo model (or “pre-statistical hypothesis’’) that 
we can construct will take into account exactly ail the factors influencing a given situation; as 
a consequence of this incomplete knowledge predictions in the social sciences are necessarily of 
a statistical nature. The last section is devoted to the problem of identification—deducing the 
value of “latent parameters” from a collection of observed data. The mathematical aspects of 
this problem were studied particularly by members of the Cowles Commission for Research in 
Economics at the University of Chicago, but Marschak shows that they can be of great im- 
portance also for the social scientist. 

In his third lecture (called “Probability and Policy”) Marschak discusses the question—and 
illustrates it by means of a simple example taken from economics—of how probabilistic models 
serve to determine and to predict actual behavior (‘‘policy” and “actions”), i.e., how they may 
be put to use in constructing “theory”’ as distinct from “mere description.” 

Marschak’s paper is written very clearly; nevertheless, the rather abstract nature of some 
sections will make it necessary to read them over carefully more than once before they can be 
thoroughly digested. 

V. “The Principal Components of Scalable Attitudes,’’ Louis Guttman. 

As pointed out by Lazarsfeld, the following two papers by Guttman are examples of 
measurement models in contrast to the models of processes previously discussed. 

For qualitative attributes which can be ranked to form a perfect scale, i.e., put into a one- 
to-one correspondence with the values of an underlying continuous variable «, Guttman had 
previously developed the mathematical theory of so-called principal components. These are 
functions of x obtained as solutions of a second-order difference equation with certain boundary 
conditions; they closely resemble, and have many qualitative properties in common with, the 
eigenfunctions of some differential equations encountered in mathematical physics. The mathe- 
matical details of the theory are not presented in this article, nor is the following practical 
question discussed: how to test whether a given set of items does form a perfect scale (at least 
approximately). The reader is referred to other publications on these subjects. Instead Gutt- 
man is now concerned primarily with the psychological meaning of these “eigenfunctions.” The 
starting point in the procedure of finding principal components is the following: to obtain a 
“new underlying metric” (i.e., a new independent variable) such that the internal consistency 
of the observed data is maximized if scale ranks are expressed in this new metric. For a perfect 
scale the first principal component fulfills this requirement. The second principal component, 
or “intensity” provides an unambiguous zero point (which we may also call “point of indiffer- 
ence’’) on the rank scale. Guttman then also proposes a psychological interpretation for the 
third and fourth principal components (called by him “closure” and ““mvolution” respective- 
ly). Clues were provided by careful analysis of experimental data, but no less so by the mathe- 
matical properties and shapes of the curves representing the principal components. This inter- 
play between psychology and mathematics is perhaps the most remarkable feature of Gutt- 
man’s paper. f : 

Whether his psychological interpretations will stand the test of further experimental veri- 
fication remains to be seen; in any event they will have provided a theory of great value for the 
stimulating and planning of further research. 

VI. “A New Approach to Factor Analysis: The Radex,” Louis Guttman. 

This model is designed explicitly to provide an explanation for the observed results of 


mental tests, although other possible applications are indicated, First Guttman discusses the 
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shortcomings of the customary procedure which consists in trying to find a (preferably small) 
number of “common factors”? to account for correlations in test scores. Then he presents his 
own model. The assumptions, stated explicitly, are that the family of tests can be sorted into 
groups of different kinds (e.g., tests on verbal ability, numerical ability, musical ability, etc.) 
and that tests of the same kind can further be ordered according to complexity (the ordering 
criterion in this case being obtained from inspection of the correlations between test results). 
Examples given in the lecture show that many series of tests can be fitted to this pattern. An 
ordered set of tests of the same kind forms what Guttman calls a perfect simplex or, rather, in 
actual practice, a quasi-simplex showing some deviations from the ideal case. A series of tests 
of the same complexity but a of different kind is arranged into a “circumflex” (i.e., a circular 
order in which no item is preferred to any other). The combination of simplex and circumflex 
then forms a structure called “radex’’ (which stands for “radial expansion of complexity’’). 

Space forbids a detailed analysis of the mathematical description which Guttman gives for 
his postulated structures (particularly the simplex and various quasi-simplexes). We can only 
mention that they are useful in predicting outcomes of further tests; in fact, Guttman writes, 
as compared to previous procedures, the radex theory offers “‘a clear path to better predictions 
with less tests.’ The relations between radex and common factor analysis are also investigated 
in detail. Unfortunately, some experimental work which was in progress at the time when 
Guttman gave his lecture and which might further have substantiated his theories could not be 
included in the book. 

Surveying both of Guttman’s contributions we agree with Lazarsfeld’s statement that they 
contain empirical findings as well as very far-reaching speculations (cf. e.g., the statement on 
p. 324 about “particle” and ‘‘wave”’ theory). It is exactly for this reason that his work deserves 
careful study and the greatest interest. 

VII. “A Conceptual Introduction to Latent Structure Analysis,’ Paul F. Lazarsfeld. 

Like Thurstone’s multiple factor hypothesis, Lazarsfeld’s latent structure analysis con- 
stitutes an attempt to explore, and to further our understanding of, the laws underlying social 
phenomena. As with Guttman’s first lecture the mathematical details of latent structure 
analysis are not given in this paper. The theory was developed in previous publications, and the 
interested reader who is not already familiar with them will certainly want to look up some of 
the references indicated by Lazarsfeld. The article under review may be called a methodological 
one: in trying to derive “latent” or “underlying”’ quantities from manifest data, some ‘‘prob- 
ability mechanisms’? must be postulated and certain further assumptions have to be intro- 
duced; these are carefully reviewed and explained by Lazarsfeld who also points out that the 
indicated procedure is not the only possible one. For simplicity this paper is limited to the case 
in which it is assumed that the latent classes can be ranked along a linear continuum. After 
having inferred a latent structure the problem arises to test whether it fits the manifest data: 
this is illustrated by an example in Section 7. Since the “raw material’’ for latent structure 
analysis consists of “itemized tests’ (i.e., in the simplest cases, collecting answers to “yes or 
no”’ questions), the relations to other techniques using the same approach are discussed in the 
final section of the paper. 


VIII. “Some Strategic Considerations in the Construction of Social Science Models,”’ 
Herbert A. Simon. 

In this paper the author, guided by his own work and experience in the field, sets forth some 
general canons of strategy which should be observed by research workers attempting to con- 
struct social science models. Simon’s program is perhaps best illustrated by quoting some of the 
questions he asks (p. 290) : “What form shall human motives take in such models, how shall the 
rational and non-rational aspects of human behavior be represented, what kind of mathe- 
matical schemes will conveniently represent the interactions of human groups, and so on?”’ 

One of Simon’s main observations is that in any particular model the specific features of 
interest arise primarily from the particular conditions under which rationality is exercised. 
Hence his third canon of strategy: “In mathematical models incorporating rational and non- 
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rational aspects of behavior, the non-rational aspects can be imbedded in the model as the 
limiting conditions that ‘bound’ the area of rational adjustment.”’ The simplest model of ra- 
tional behavior is a static one, based on the individual’s maximizing a certain “criterion” (the 
latter being a quantity incorporated in the model). But as our considerations are made more 
realistic, taking into account the dynamics of the individual’s adjustment and the fact that one 
person is confronted by the varying and, to some extent, unpredictable reactions of others, we 
have a shift from the model of optimization toward what Simon calls the model of adaptive 
behavior. 

It may be permitted to point out here what looks like a minor contradiction in Simon’s 
paper. He writes at the beginning (p. 388) : “Mathematical social science is, first and foremost, 
social science. If it is bad social science (i.e., empirically false) the fact that it is good mathe- 
matics (i.e., logically consistent) should provide little comfort.’’ But later (p. 409) when he 
gives a mathematical description of Homans’ “Human Group” he states: ‘“‘. . . the reader can 
judge for himself how well we have succeeded in capturing the essential features of Homans’ 
system in our equations. In any event it is unnecessary to concern ourselves here with the 
exactness of the representation or the empirical correctness of his postulates.’’ In other words, 
we are asked to assume that Homans’ theory is “empirically correct,” and, even granting that, 
no effort is undertaken to make sure that the system of equations derived by Simon represents 
‘good social science.” It is, however, true that the model’s validity can be verified a posteriori, 
by showing that some of the conclusions drawn from it agree with empirical results. 

After this brief survey of the eight articles, we should direct the reader’s attention to the 
various interrelations and cross references among them. Some are quite obvious, like, e.g., the 
connection between Marschak’s problem of identification and Lazarsfeld’s latent structure 
analysis. Others will emerge as one rereads the papers. For instance, one might try to apply 
Coleman’s considerations to the models of Anderson and Simon, or to follow up some of the 
hints contained in Guttman’s lectures. Further relations are outlined in Lazarsfeld’s introduc- 
tion; they are, of course, not revealed by a superficial glancing at the book; attentive reading 
will be found very profitable. 

It should be made clear that most of the material presented is not really new, the most 
notable exception being Guttman’s psychological interpretation of the third and fourth princi- 
pal components, which is given here in detail for the first time. The rest had either been pub- 
lished elsewhere (and often in greater detail than would be possible in a short lecture), or else it 
constitutes an able presentation and analysis of already existing and well-known material. How- 
ever, it is certainly quite useful to have various lines of scientific inquiry collected in one vol- 
ume; this book will prove a good source of information and ideas to anyone interested in theo- 
retical social science. 

Unfortunately, there are, in this book, far too many misprints which, in places, even tend to 
make the meaning of the text unintelligible. It is the reviewer’s understanding that this hap- 
pened because the individual contributors were not given the opportunity to read the proofs 
of their contributions. We certainly hope that in later editions of the book most of the misprints 


will be eliminated. 
ERNESTO TRUCCO 
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